N. gono Cost of Resistance

set.seed(1) 
library(splines)
library(ggplot2)
library(ggtree)
#> ggtree v3.6.2 For help: https://yulab-smu.top/treedata-book/
#> 
#> If you use the ggtree package suite in published research, please cite
#> the appropriate paper(s):
#> 
#> Guangchuang Yu, David Smith, Huachen Zhu, Yi Guan, Tommy Tsan-Yuk Lam.
#> ggtree: an R package for visualization and annotation of phylogenetic
#> trees with their covariates and other associated data. Methods in
#> Ecology and Evolution. 2017, 8(1):28-36. doi:10.1111/2041-210X.12628
#> 
#> Guangchuang Yu.  Data Integration, Manipulation and Visualization of
#> Phylogenetic Trees (1st edition). Chapman and Hall/CRC. 2022,
#> doi:10.1201/9781003279242
#> 
#> Guangchuang Yu, Tommy Tsan-Yuk Lam, Huachen Zhu, Yi Guan. Two methods
#> for mapping and visualizing associated data on phylogeny using ggtree.
#> Molecular Biology and Evolution. 2018, 35(12):3041-3043.
#> doi:10.1093/molbev/msy194
library(reshape2)
library(ResistPhy)
library(ape)
#> 
#> Attaching package: 'ape'
#> The following object is masked from 'package:ggtree':
#> 
#>     rotate
library(posterior)
#> This is posterior version 1.4.0
#> 
#> Attaching package: 'posterior'
#> The following objects are masked from 'package:stats':
#> 
#>     mad, sd, var
#> The following objects are masked from 'package:base':
#> 
#>     %in%, match
library(cmdstanr)
#> This is cmdstanr version 0.5.3
#> - CmdStanR documentation and vignettes: mc-stan.org/cmdstanr
#> - CmdStan path: /home/david/.cmdstan/cmdstan-2.31.0
#> - CmdStan version: 2.31.0
library(patchwork)
library(viridis)
#> Loading required package: viridisLite
library(RColorBrewer)
library(latex2exp)

run_mcmc <- F

Load tree from Yonathan Grad’s paper, do some boring bioinformatics stuff

Honestly wish my computer could just read my mind and do what I want, automatically

tr <- read.tree(system.file("extdata", "grad2016.nwk", package="ResistPhy", mustWork=T))
tr <- makeNodeLabel(tr)
meta <- read.table(system.file("extdata", "grad2016.tab", package="ResistPhy", mustWork=T),
            sep='\t',comment.char='',header=T,as.is=T)
rownames(meta) <- meta$ID
cutoffs <- read.csv(system.file("extdata", "cutoffs.csv", package="ResistPhy", mustWork=T))

Dichotomise Resistance data

rdf <- data.frame(ID=meta$ID)
rownames(rdf) <- rownames(meta)

abx_names <- data.frame(cutoff_name = c("CFX", "CIP", "CRO", "AZI", "PEN", "TET", "SPC"),
    meta_name=c("CFX", "CIP", "CRO", "AZI", "PEN", "TET", "SPC")
)

for (s in c(1:nrow(abx_names))) {
    rdf[[abx_names$cutoff_name[s]]] <- (meta[[abx_names$meta_name[s]]] > cutoffs[[abx_names$cutoff_name[s]]][2])
}

rdf<-within(rdf, rm(ID))

Select the three lineages of interest

CL_1 <- extract.clade(tr, "Node129")$tip.label
CL_2 <- extract.clade(tr, "Node580")$tip.label
CL_3 <- extract.clade(tr, "Node1038")$tip.label

tr_r <- keep.tip(tr, c(CL_1,CL_2,CL_3))

#Remove CFX resistant subclade
tr_r <- drop.tip(tr_r, extract.clade(tr_r, "Node732")$tip.label)

tr_r <- ladderize(tr_r)

p <- ggtree(tr_r) +  geom_nodelab(geom='text',size=4) +
    theme_tree2()
gheatmap(p, rdf, offset=8, width=0.6, 
        colnames=FALSE, legend_title="Resistance")+
        scale_x_ggtree()
#> Scale for y is already present.
#> Adding another scale for y, which will replace the existing scale.

We are interested in the two lineages within the subtree with root at Node125

cl_dat <- data.frame(node=nodeid(tr_r, c("Node580", "Node129", "Node1038")), name=c("FQ Susceptible", "FQ Resistant 1","FQ Resistant 2"))
mrsd_full <- max(meta[tr_r$tip.label, "Year"])

p <- ggtree(tr_r, mrsd=paste0(mrsd_full,"-01-01")) + 
    geom_cladelab(
        data = cl_dat,
        mapping = aes(
            node = node, 
            label = name, 
            color = name
        ),
        fontsize = 3,
        angle=90,
        align = TRUE,
        vjust = 0.0,
        hjust = "center",
        offset= 10,
        lwd=4.0,
        show.legend = FALSE
    ) + 
    scale_color_brewer(palette="Dark2") +
    theme_tree2()
h1 <- gheatmap(p, rdf[,c("CIP", "AZI", "PEN", "TET")], width=0.2, colnames=FALSE, legend_title="Resistance")+
        scale_fill_manual(values=c("blue", "red")) + 
        guides(color = "none", fill=guide_legend("MIC >= Cutoff")) +
        scale_x_ggtree() +
        theme(legend.position="bottom",
        axis.text.x = element_text(angle = 45, vjust = 0.5, hjust=1))
#> Scale for y is already present.
#> Adding another scale for y, which will replace the existing scale.
#> Scale for fill is already present.
#> Adding another scale for fill, which will replace the existing scale.

h1

pdf("Figures/gono_tree_90.pdf",8,8)
plot(h1)
dev.off()
#> png 
#>   2

Node580 corresponds to FQ susceptible, Node129 is FQ resistant Extract the two lineages, filter out tips with missing FQ profiles,or that are FQ resistant and in the susceptible clade and vice versa Further remove all tips that are CIP resistant (mostly a small subclade in the susceptible tree)

sus_tr <- extract.clade(tr_r, "Node580")
res_tr_1 <- extract.clade(tr_r, "Node129")
res_tr_2 <- extract.clade(tr_r, "Node1038")

rdf_sus <- rdf[sus_tr$tip.label,]
rdf_res_1 <- rdf[res_tr_1$tip.label,]
rdf_res_2 <- rdf[res_tr_2$tip.label,]

tip_times_sus <- meta[rownames(rdf_sus),]$Year
tip_times_res_1 <- meta[rownames(rdf_res_1),]$Year
tip_times_res_2 <- meta[rownames(rdf_res_2),]$Year

mrsd_sus <- max(tip_times_sus)
mrsd_res_1 <- max(tip_times_res_1)
mrsd_res_2 <- max(tip_times_res_2)

Plot the filtered subtrees

p1 <- ggtree(sus_tr, mrsd=paste0(mrsd_sus,"-01-01")) +
    geom_tiplab(size=1, align=TRUE, linesize=.5)+
    scale_x_ggtree() + 
    theme_tree2()
gheatmap(p1, rdf, offset=8, width=0.6, 
    colnames=FALSE, legend_title="Resistance")+
    scale_fill_manual(values=c("blue", "red"))+
    labs(title="Susceptible Tree") +
    guides(color = "none", fill=guide_legend("MIC >= Cutoff")) +
    scale_x_ggtree()
#> Scale for y is already present.
#> Adding another scale for y, which will replace the existing scale.
#> Scale for fill is already present.
#> Adding another scale for fill, which will replace the existing scale.
#> Scale for x is already present.
#> Adding another scale for x, which will replace the existing scale.

p2 <-ggtree(res_tr_1, mrsd=paste0(mrsd_res_1,"-01-01"))+theme_tree2() +
    geom_tiplab(size=1, align=TRUE, linesize=.5)+
    scale_x_ggtree() + 
    theme_tree2()
gheatmap(p2, rdf, offset=8, width=0.6, 
        colnames=FALSE, legend_title="Resistance")+
        scale_fill_manual(values=c("blue", "red"))+
        labs(title="Resistant Tree") +
        guides(color = "none", fill=guide_legend("MIC >= Cutoff")) +
        scale_x_ggtree()
#> Scale for y is already present.
#> Adding another scale for y, which will replace the existing scale.
#> Scale for fill is already present.
#> Adding another scale for fill, which will replace the existing scale.
#> Scale for x is already present.
#> Adding another scale for x, which will replace the existing scale.

p2 <-ggtree(res_tr_2, mrsd=paste0(mrsd_res_2,"-01-01"))+theme_tree2() +
    geom_tiplab(size=1, align=TRUE, linesize=.5)+
    scale_x_ggtree() + 
    theme_tree2()
gheatmap(p2, rdf, offset=8, width=0.6, 
        colnames=FALSE, legend_title="Resistance")+
        scale_fill_manual(values=c("blue", "red"))+
        labs(title="Resistant Tree") +
        guides(color = "none", fill=guide_legend("MIC >= Cutoff")) +
        scale_x_ggtree()
#> Scale for y is already present.
#> Adding another scale for y, which will replace the existing scale.
#> Scale for fill is already present.
#> Adding another scale for fill, which will replace the existing scale.
#> Scale for x is already present.
#> Adding another scale for x, which will replace the existing scale.

Set the analysis end date, move tips into the middle of the year as opposed to the beginning by offsetting most recent sampling date by +0.5

t_end <- max(mrsd_res_1, mrsd_res_2, mrsd_sus)+1 
mr_sus <- t_end-mrsd_sus-0.5
mr_res_1 <- t_end-mrsd_res_1-0.5
mr_res_2 <- t_end-mrsd_res_2-0.5

print(t_end)
#> [1] 2014

Load and preprocess usage data

pl_begin <- 1988
# AZI, PEN, TET, CEF, FQ
dat <- system.file("extdata", "usage-clean.csv", package = "ResistPhy", mustWork = TRUE)
abx_raw <- as.matrix(read.csv(dat, header=F))
abx_df <- as.data.frame(t(abx_raw))
colnames(abx_df) <- c("year",abx_df[1,2:14])
abx_df <- abx_df[-1,]
rownames(abx_df) <- as.integer(abx_df[,1])
abx_df$year <- as.integer(abx_df$year)
abx_df$FQ <- as.numeric(abx_df[, "Ciprofloxacin"]) + as.numeric(abx_df[, "Ofloxacin"])
abx_df$CPH <- as.numeric(abx_df[, "Cefixime"]) + as.numeric(abx_df[, "Ceftriaxone 250 mg"]) + as.numeric(abx_df[, "Ceftriaxone 125 mg"]) +as.numeric(abx_df[, "Other Cephalo."]) 

abx_df$Other <- 100-as.numeric(abx_df[, "FQ"]) - 
                as.numeric(abx_df[, "CPH"]) -
                as.numeric(abx_df[, "Penicillins"]) -
                as.numeric(abx_df[, "Azithromycin 2gm"]) -
                as.numeric(abx_df[, "Tetracyclines"])

usage_f_fq <- yearly_usg_stepfunc(as.numeric(abx_df[sapply(c(pl_begin:t_end), paste0),"FQ"]), c(pl_begin:t_end))
usage_f_pen <- yearly_usg_stepfunc(as.numeric(abx_df[sapply(c(pl_begin:t_end), paste0),"Penicillins"]), c(pl_begin:t_end))
usage_f_tet <- yearly_usg_stepfunc(as.numeric(abx_df[sapply(c(pl_begin:t_end), paste0),"Tetracyclines"]), c(pl_begin:t_end))
usage_f_azi <- yearly_usg_stepfunc(as.numeric(abx_df[sapply(c(pl_begin:t_end), paste0),"Azithromycin 2gm"]), c(pl_begin:t_end))
usage_f_CPH <- yearly_usg_stepfunc(as.numeric(abx_df[sapply(c(pl_begin:t_end), paste0),"CPH"]), c(pl_begin:t_end))
usage_f_other <- yearly_usg_stepfunc(as.numeric(abx_df[sapply(c(pl_begin:t_end), paste0),"Other"]), c(pl_begin:t_end))

time_grid <- seq(from=pl_begin, to=t_end, length.out=4000)
fq_df <- data.frame(year=rep(time_grid,6), 
                    usage=c(sapply(time_grid,usage_f_fq),
                            sapply(time_grid,usage_f_pen),
                            sapply(time_grid,usage_f_tet),
                            sapply(time_grid,usage_f_azi),
                            sapply(time_grid,usage_f_CPH),
                            sapply(time_grid,usage_f_other)),
                    abx=c(rep("Fluoroquinolones",length(time_grid)), 
                          rep("Penicillins", length(time_grid)), 
                          rep("Tetracyclines",length(time_grid)),
                          rep("Azithromycin", length(time_grid)),
                          rep("Cephalosporins", length(time_grid)),
                          rep("Other", length(time_grid))))

Plot usage data, plot individual areas for AMRs that the lineages carry resistance against

p <- ggplot(fq_df, aes(x=year, y=usage, fill=factor(abx, levels=rev(c("Fluoroquinolones", "Penicillins","Tetracyclines", "Cephalosporins", "Azithromycin", "Other"))))) +
    geom_area() +
    scale_x_continuous(breaks = seq(pl_begin, t_end, by =1)) + 
    labs(x="Year", y="Average Usage as Proportion of Primary Treatment") +
    guides(fill = guide_legend(title="Antimicrobial",reverse=F)) +
    scale_fill_viridis(discrete=T) +
    scale_y_continuous(labels = scales::percent_format(scale = 1), breaks=seq(from=0,to=100, by=10)) +
    geom_vline(xintercept=1995, linetype="longdash", color="gray70",alpha=.8)+
    geom_text(aes(x=1995, label="\nAnalysis Start Date", y=40), size=rel(4.0), colour="gray70", alpha=.8, angle=90)+
    theme_minimal() +
    theme(axis.text.x=element_text(size=rel(0.7), angle = 45, vjust=0.5, hjust=0.5),
        axis.text.y=element_text(size=rel(0.7), hjust=1),
        axis.title.y=element_text(size=rel(0.6)),
        axis.title.x=element_text(size=rel(0.6)),
        aspect.ratio=1,
        panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(),
        axis.line = element_line(size=rel(0.2), colour = "grey80"),
        legend.position="right",
        legend.justification="left",
        legend.margin=margin(0,0,0,0),
        legend.box.margin=margin(0,0,0,-10))
#> Warning: The `size` argument of `element_line()` is deprecated as of ggplot2 3.4.0.
#> ℹ Please use the `linewidth` argument instead.

pl1 <- h1 + annotation_custom(
    grob = ggplotGrob(p),
    xmin = 1795,
    xmax = 1995,
    ymin = 190,
    ymax = 410)

pdf("Figures/gono_data.pdf",6,6)
pl1
dev.off()
#> png 
#>   2

We can see that PEN + TET was used in significant quantities before ~1995. The FQ resistant lineage carries resistance against these as well, but susceptible clade does not. Set the analysis start date to 1995. Subset FQ usage data for the given timespan.

t_begin <- 1995
usage_df <- data.frame(year=c(t_begin:t_end), usage=as.numeric(abx_df[sapply(c(t_begin:t_end), paste0),"FQ"]))
usage_df$usage <- usage_df$usage/100

Perform Inference

• Time scale is in years
• Infectious period of gono is ~12 days

gamma_ansatz <- 1/90.0
time_scale <- 365.0
t0 <- t_begin #0
tmax <- t_end #t_end - t_begin
usage_df$time #<- usage_df$time - t_begin
#> NULL
stopifnot(all(usage_df$time<=tmax))
stopifnot(all(usage_df$time>=t0))

Sample

out <- infer_costs2(list(sus_tr, res_tr_1, res_tr_2), 
                    list(mr_sus,mr_res_1, mr_res_2),
                    usage_df$usage,
                    usage_df$year,
                    t0,
                    tmax,
                    gamma_ansatz,
                    time_scale,
                    n_iter=2000, 
                    n_warmup=2000,
                    model="nodecay",
                    K=60,
                    L=6.5,
                    seed=1345678,
                    gamma_log_sd = 0.1,
                    stan_control=list(adapt_delta=.99,
                        max_treedepth=13,
                        parallel_chains=4,
                        chains=4,
                        refresh=1000
                    ))
saveRDS(out, "~/gono_90s.rds")

Load existing analysis result

out <- readRDS("~/gono_90s.rds")

Check convergence

print(out$converged)
#> [1] TRUE

Inferred dynamics & Rt

source("Figures/figure3.R")
pdf("Figures/figure3_90.pdf",6,6)
plot(fig3)
dev.off()
#> png 
#>   2
fig3

Pairs Plot

fig4 <- plot_qpairs(out)
pdf("Figures/figure4_90.pdf",6,6)
plot(fig4)
dev.off()
#> png 
#>   2
fig4

RR map

source("Figures/figure5.R")
#> Warning: Removed 7200 rows containing non-finite values
#> (`stat_contour_filled()`).
#> Warning: Removed 7800 rows containing non-finite values
#> (`stat_contour_filled()`).
#> Warning: Removed 7200 rows containing non-finite values
#> (`stat_contour_filled()`).
#> Warning: Removed 7800 rows containing non-finite values
#> (`stat_contour_filled()`).
#> Scale for y is already present.
#> Adding another scale for y, which will replace the existing scale.
#> Scale for y is already present.
#> Adding another scale for y, which will replace the existing scale.
#> Warning: Removed 7200 rows containing non-finite values (`stat_contour_filled()`).
#> Removed 7800 rows containing non-finite values (`stat_contour_filled()`).
#> Warning: Removed 7200 rows containing non-finite values
#> (`stat_contour_filled()`).
#> Warning: Removed 7800 rows containing non-finite values
#> (`stat_contour_filled()`).
#> Scale for x is already present.
#> Adding another scale for x, which will replace the existing scale.
#> Scale for x is already present.
#> Adding another scale for x, which will replace the existing scale.
pdf("Figures/figure5_90.pdf",6,6)
plot(fig5)
#> Warning: Removed 7200 rows containing non-finite values
#> (`stat_contour_filled()`).
#> Warning: Removed 7200 rows containing non-finite values
#> (`stat_contour_filled()`).
dev.off()
#> png 
#>   2
fig5
#> Warning: Removed 7200 rows containing non-finite values (`stat_contour_filled()`).
#> Removed 7200 rows containing non-finite values (`stat_contour_filled()`).

Traces & ppcheck

s4 <- plot_traces(out)
pdf("Figures/s4_90.pdf",12,12)
plot(s4)
dev.off()
#> png 
#>   2
s4

s5a <- plot_ppcheck_At(out, 1)
#>  [1]  0.5  1.5  2.5  3.5  4.5  5.5  6.5  7.5  8.5  9.5 10.5 11.5 12.5 13.5
#>  [1] 35 26 13 23 15  1 10  5 18 16  8  6  5  9
#> [1] 2014
#>    [1] 1995.000 1995.009 1995.019 1995.028 1995.038 1995.047 1995.056 1995.066
#>    [9] 1995.075 1995.084 1995.094 1995.103 1995.113 1995.122 1995.131 1995.141
#>   [17] 1995.150 1995.159 1995.169 1995.178 1995.188 1995.197 1995.206 1995.216
#>   [25] 1995.225 1995.235 1995.244 1995.253 1995.263 1995.272 1995.281 1995.291
#>   [33] 1995.300 1995.310 1995.319 1995.328 1995.338 1995.347 1995.356 1995.366
#>   [41] 1995.375 1995.385 1995.394 1995.403 1995.413 1995.422 1995.432 1995.441
#>   [49] 1995.450 1995.460 1995.469 1995.478 1995.488 1995.497 1995.507 1995.516
#>   [57] 1995.525 1995.535 1995.544 1995.553 1995.563 1995.572 1995.582 1995.591
#>   [65] 1995.600 1995.610 1995.619 1995.628 1995.638 1995.647 1995.656 1995.666
#>   [73] 1995.675 1995.685 1995.694 1995.703 1995.713 1995.722 1995.731 1995.741
#>   [81] 1995.750 1995.760 1995.769 1995.778 1995.788 1995.797 1995.807 1995.816
#>   [89] 1995.825 1995.835 1995.844 1995.854 1995.863 1995.873 1995.882 1995.891
#>   [97] 1995.901 1995.910 1995.920 1995.929 1995.938 1995.948 1995.957 1995.967
#>  [105] 1995.976 1995.986 1995.995 1996.004 1996.014 1996.023 1996.033 1996.042
#>  [113] 1996.051 1996.061 1996.070 1996.080 1996.089 1996.098 1996.108 1996.117
#>  [121] 1996.127 1996.136 1996.146 1996.155 1996.164 1996.174 1996.183 1996.193
#>  [129] 1996.202 1996.211 1996.221 1996.230 1996.240 1996.249 1996.259 1996.268
#>  [137] 1996.277 1996.287 1996.296 1996.306 1996.315 1996.324 1996.334 1996.343
#>  [145] 1996.353 1996.362 1996.371 1996.381 1996.390 1996.399 1996.409 1996.418
#>  [153] 1996.428 1996.437 1996.446 1996.456 1996.465 1996.475 1996.484 1996.493
#>  [161] 1996.503 1996.512 1996.522 1996.531 1996.540 1996.550 1996.559 1996.568
#>  [169] 1996.578 1996.587 1996.597 1996.606 1996.615 1996.625 1996.634 1996.644
#>  [177] 1996.653 1996.662 1996.672 1996.681 1996.691 1996.700 1996.709 1996.719
#>  [185] 1996.728 1996.737 1996.747 1996.756 1996.766 1996.775 1996.784 1996.794
#>  [193] 1996.803 1996.813 1996.822 1996.831 1996.841 1996.850 1996.860 1996.869
#>  [201] 1996.878 1996.888 1996.897 1996.907 1996.916 1996.926 1996.935 1996.944
#>  [209] 1996.954 1996.963 1996.973 1996.982 1996.992 1997.001 1997.010 1997.020
#>  [217] 1997.029 1997.039 1997.048 1997.058 1997.067 1997.077 1997.086 1997.095
#>  [225] 1997.105 1997.114 1997.124 1997.133 1997.143 1997.152 1997.161 1997.171
#>  [233] 1997.180 1997.190 1997.199 1997.209 1997.218 1997.228 1997.237 1997.246
#>  [241] 1997.256 1997.265 1997.275 1997.284 1997.294 1997.303 1997.313 1997.322
#>  [249] 1997.331 1997.341 1997.350 1997.360 1997.369 1997.379 1997.388 1997.397
#>  [257] 1997.407 1997.416 1997.426 1997.435 1997.445 1997.454 1997.464 1997.473
#>  [265] 1997.482 1997.492 1997.501 1997.511 1997.520 1997.530 1997.539 1997.548
#>  [273] 1997.558 1997.567 1997.577 1997.586 1997.596 1997.605 1997.615 1997.624
#>  [281] 1997.633 1997.643 1997.652 1997.662 1997.671 1997.681 1997.690 1997.699
#>  [289] 1997.709 1997.718 1997.728 1997.737 1997.747 1997.756 1997.766 1997.775
#>  [297] 1997.784 1997.794 1997.803 1997.813 1997.822 1997.832 1997.841 1997.850
#>  [305] 1997.860 1997.869 1997.879 1997.888 1997.898 1997.907 1997.917 1997.926
#>  [313] 1997.935 1997.945 1997.953 1997.961 1997.969 1997.978 1997.987 1997.996
#>  [321] 1998.005 1998.015 1998.024 1998.033 1998.042 1998.051 1998.061 1998.070
#>  [329] 1998.079 1998.088 1998.097 1998.107 1998.116 1998.125 1998.134 1998.143
#>  [337] 1998.153 1998.162 1998.171 1998.180 1998.189 1998.199 1998.208 1998.217
#>  [345] 1998.226 1998.235 1998.244 1998.253 1998.262 1998.272 1998.281 1998.290
#>  [353] 1998.299 1998.308 1998.318 1998.327 1998.336 1998.345 1998.355 1998.364
#>  [361] 1998.373 1998.382 1998.392 1998.401 1998.410 1998.419 1998.428 1998.438
#>  [369] 1998.447 1998.456 1998.465 1998.475 1998.484 1998.493 1998.502 1998.511
#>  [377] 1998.521 1998.530 1998.539 1998.548 1998.558 1998.567 1998.576 1998.585
#>  [385] 1998.587 1998.596 1998.606 1998.615 1998.624 1998.633 1998.643 1998.652
#>  [393] 1998.661 1998.671 1998.680 1998.689 1998.698 1998.708 1998.717 1998.726
#>  [401] 1998.735 1998.745 1998.754 1998.763 1998.773 1998.782 1998.791 1998.800
#>  [409] 1998.810 1998.819 1998.828 1998.838 1998.847 1998.856 1998.865 1998.875
#>  [417] 1998.884 1998.893 1998.902 1998.912 1998.921 1998.930 1998.940 1998.949
#>  [425] 1998.958 1998.967 1998.976 1998.986 1998.995 1999.004 1999.013 1999.020
#>  [433] 1999.027 1999.036 1999.046 1999.055 1999.065 1999.074 1999.084 1999.093
#>  [441] 1999.103 1999.112 1999.122 1999.131 1999.141 1999.150 1999.160 1999.169
#>  [449] 1999.179 1999.188 1999.198 1999.207 1999.217 1999.226 1999.236 1999.245
#>  [457] 1999.255 1999.264 1999.274 1999.283 1999.292 1999.302 1999.311 1999.321
#>  [465] 1999.330 1999.340 1999.349 1999.358 1999.367 1999.376 1999.385 1999.394
#>  [473] 1999.403 1999.412 1999.421 1999.430 1999.438 1999.447 1999.456 1999.465
#>  [481] 1999.474 1999.482 1999.491 1999.500 1999.509 1999.519 1999.528 1999.537
#>  [489] 1999.546 1999.556 1999.565 1999.574 1999.584 1999.593 1999.602 1999.611
#>  [497] 1999.620 1999.629 1999.639 1999.648 1999.657 1999.666 1999.675 1999.684
#>  [505] 1999.694 1999.703 1999.712 1999.721 1999.731 1999.740 1999.749 1999.758
#>  [513] 1999.768 1999.777 1999.786 1999.796 1999.805 1999.814 1999.823 1999.833
#>  [521] 1999.842 1999.851 1999.860 1999.870 1999.879 1999.888 1999.898 1999.907
#>  [529] 1999.916 1999.925 1999.935 1999.944 1999.953 1999.963 1999.972 1999.981
#>  [537] 1999.990 2000.000 2000.009 2000.018 2000.027 2000.037 2000.046 2000.055
#>  [545] 2000.065 2000.074 2000.083 2000.093 2000.102 2000.111 2000.120 2000.130
#>  [553] 2000.139 2000.148 2000.158 2000.167 2000.176 2000.186 2000.195 2000.204
#>  [561] 2000.213 2000.223 2000.232 2000.241 2000.250 2000.260 2000.269 2000.278
#>  [569] 2000.287 2000.297 2000.306 2000.315 2000.324 2000.334 2000.343 2000.352
#>  [577] 2000.361 2000.371 2000.380 2000.389 2000.398 2000.408 2000.417 2000.426
#>  [585] 2000.435 2000.445 2000.454 2000.463 2000.472 2000.482 2000.491 2000.500
#>  [593] 2000.508 2000.515 2000.523 2000.531 2000.539 2000.548 2000.557 2000.566
#>  [601] 2000.573 2000.580 2000.588 2000.595 2000.605 2000.614 2000.623 2000.632
#>  [609] 2000.641 2000.651 2000.660 2000.669 2000.678 2000.687 2000.697 2000.706
#>  [617] 2000.715 2000.724 2000.733 2000.743 2000.752 2000.761 2000.770 2000.779
#>  [625] 2000.789 2000.798 2000.807 2000.816 2000.825 2000.835 2000.844 2000.853
#>  [633] 2000.863 2000.872 2000.882 2000.891 2000.900 2000.910 2000.919 2000.928
#>  [641] 2000.938 2000.947 2000.957 2000.966 2000.975 2000.985 2000.994 2001.004
#>  [649] 2001.013 2001.022 2001.032 2001.041 2001.050 2001.060 2001.069 2001.079
#>  [657] 2001.087 2001.096 2001.104 2001.113 2001.121 2001.130 2001.138 2001.148
#>  [665] 2001.157 2001.167 2001.176 2001.185 2001.195 2001.204 2001.214 2001.223
#>  [673] 2001.232 2001.241 2001.251 2001.260 2001.269 2001.278 2001.287 2001.297
#>  [681] 2001.306 2001.315 2001.324 2001.334 2001.343 2001.352 2001.361 2001.371
#>  [689] 2001.380 2001.389 2001.398 2001.408 2001.417 2001.426 2001.435 2001.445
#>  [697] 2001.454 2001.463 2001.472 2001.482 2001.491 2001.500 2001.506 2001.515
#>  [705] 2001.525 2001.534 2001.543 2001.553 2001.562 2001.572 2001.581 2001.590
#>  [713] 2001.599 2001.608 2001.618 2001.627 2001.636 2001.645 2001.654 2001.664
#>  [721] 2001.673 2001.682 2001.691 2001.700 2001.710 2001.719 2001.727 2001.736
#>  [729] 2001.745 2001.753 2001.762 2001.769 2001.777 2001.784 2001.791 2001.800
#>  [737] 2001.810 2001.819 2001.828 2001.837 2001.846 2001.855 2001.865 2001.874
#>  [745] 2001.883 2001.892 2001.901 2001.911 2001.920 2001.929 2001.938 2001.947
#>  [753] 2001.957 2001.966 2001.975 2001.984 2001.993 2002.002 2002.012 2002.021
#>  [761] 2002.030 2002.040 2002.049 2002.058 2002.068 2002.077 2002.086 2002.096
#>  [769] 2002.105 2002.115 2002.124 2002.133 2002.143 2002.152 2002.161 2002.171
#>  [777] 2002.180 2002.189 2002.198 2002.208 2002.217 2002.226 2002.235 2002.244
#>  [785] 2002.254 2002.263 2002.272 2002.281 2002.290 2002.300 2002.309 2002.318
#>  [793] 2002.327 2002.336 2002.346 2002.355 2002.364 2002.373 2002.382 2002.392
#>  [801] 2002.400 2002.408 2002.417 2002.425 2002.433 2002.442 2002.450 2002.458
#>  [809] 2002.462 2002.469 2002.475 2002.482 2002.491 2002.500 2002.509 2002.519
#>  [817] 2002.528 2002.538 2002.547 2002.556 2002.564 2002.573 2002.581 2002.590
#>  [825] 2002.599 2002.609 2002.618 2002.627 2002.637 2002.646 2002.655 2002.665
#>  [833] 2002.674 2002.683 2002.693 2002.702 2002.712 2002.721 2002.731 2002.740
#>  [841] 2002.750 2002.759 2002.769 2002.778 2002.788 2002.797 2002.807 2002.816
#>  [849] 2002.826 2002.835 2002.845 2002.854 2002.864 2002.873 2002.883 2002.892
#>  [857] 2002.902 2002.911 2002.920 2002.930 2002.939 2002.949 2002.957 2002.965
#>  [865] 2002.975 2002.984 2002.993 2003.003 2003.012 2003.021 2003.031 2003.040
#>  [873] 2003.049 2003.059 2003.068 2003.077 2003.087 2003.095 2003.104 2003.112
#>  [881] 2003.120 2003.130 2003.139 2003.148 2003.157 2003.166 2003.175 2003.184
#>  [889] 2003.193 2003.203 2003.212 2003.221 2003.230 2003.239 2003.248 2003.254
#>  [897] 2003.260 2003.266 2003.271 2003.272 2003.281 2003.290 2003.299 2003.308
#>  [905] 2003.317 2003.326 2003.335 2003.345 2003.354 2003.363 2003.372 2003.381
#>  [913] 2003.390 2003.399 2003.409 2003.418 2003.427 2003.436 2003.445 2003.454
#>  [921] 2003.463 2003.473 2003.482 2003.491 2003.500 2003.509 2003.516 2003.524
#>  [929] 2003.533 2003.543 2003.552 2003.561 2003.570 2003.579 2003.589 2003.598
#>  [937] 2003.607 2003.616 2003.625 2003.634 2003.643 2003.652 2003.661 2003.670
#>  [945] 2003.679 2003.688 2003.697 2003.705 2003.712 2003.720 2003.727 2003.735
#>  [953] 2003.743 2003.750 2003.758 2003.765 2003.771 2003.778 2003.787 2003.797
#>  [961] 2003.806 2003.815 2003.824 2003.833 2003.843 2003.852 2003.861 2003.870
#>  [969] 2003.879 2003.889 2003.898 2003.907 2003.916 2003.925 2003.935 2003.944
#>  [977] 2003.951 2003.959 2003.967 2003.976 2003.985 2003.995 2004.004 2004.014
#>  [985] 2004.023 2004.032 2004.041 2004.050 2004.059 2004.067 2004.075 2004.083
#>  [993] 2004.092 2004.100 2004.108 2004.116 2004.123 2004.130 2004.137 2004.143
#> [1001] 2004.152 2004.161 2004.169 2004.178 2004.187 2004.195 2004.204 2004.213
#> [1009] 2004.221 2004.230 2004.237 2004.244 2004.251 2004.260 2004.269 2004.278
#> [1017] 2004.287 2004.297 2004.306 2004.315 2004.324 2004.333 2004.342 2004.351
#> [1025] 2004.360 2004.369 2004.379 2004.388 2004.397 2004.406 2004.416 2004.425
#> [1033] 2004.434 2004.444 2004.453 2004.463 2004.472 2004.481 2004.491 2004.500
#> [1041] 2004.509 2004.519 2004.528 2004.538 2004.547 2004.557 2004.566 2004.576
#> [1049] 2004.585 2004.595 2004.604 2004.613 2004.622 2004.631 2004.640 2004.649
#> [1057] 2004.658 2004.666 2004.673 2004.681 2004.688 2004.696 2004.704 2004.712
#> [1065] 2004.720 2004.728 2004.735 2004.743 2004.752 2004.761 2004.771 2004.780
#> [1073] 2004.790 2004.799 2004.808 2004.818 2004.827 2004.837 2004.846 2004.855
#> [1081] 2004.865 2004.874 2004.884 2004.893 2004.903 2004.912 2004.921 2004.931
#> [1089] 2004.940 2004.950 2004.959 2004.968 2004.978 2004.987 2004.997 2005.006
#> [1097] 2005.016 2005.025 2005.034 2005.044 2005.053 2005.063 2005.072 2005.081
#> [1105] 2005.091 2005.100 2005.110 2005.119 2005.129 2005.138 2005.147 2005.157
#> [1113] 2005.166 2005.176 2005.185 2005.194 2005.204 2005.213 2005.223 2005.229
#> [1121] 2005.236 2005.242 2005.249 2005.256 2005.265 2005.274 2005.284 2005.293
#> [1129] 2005.302 2005.311 2005.320 2005.330 2005.339 2005.348 2005.357 2005.366
#> [1137] 2005.376 2005.385 2005.394 2005.403 2005.412 2005.422 2005.431 2005.440
#> [1145] 2005.449 2005.458 2005.468 2005.477 2005.486 2005.493 2005.500 2005.509
#> [1153] 2005.519 2005.528 2005.537 2005.546 2005.556 2005.565 2005.574 2005.583
#> [1161] 2005.593 2005.602 2005.611 2005.620 2005.630 2005.639 2005.648 2005.657
#> [1169] 2005.667 2005.676 2005.685 2005.695 2005.704 2005.713 2005.723 2005.732
#> [1177] 2005.741 2005.751 2005.760 2005.769 2005.779 2005.788 2005.797 2005.806
#> [1185] 2005.816 2005.825 2005.834 2005.844 2005.853 2005.862 2005.867 2005.876
#> [1193] 2005.885 2005.894 2005.902 2005.911 2005.920 2005.929 2005.938 2005.946
#> [1201] 2005.955 2005.964 2005.973 2005.982 2005.991 2006.000 2006.010 2006.019
#> [1209] 2006.027 2006.035 2006.044 2006.053 2006.063 2006.072 2006.081 2006.090
#> [1217] 2006.099 2006.108 2006.117 2006.126 2006.135 2006.144 2006.154 2006.163
#> [1225] 2006.172 2006.181 2006.190 2006.199 2006.208 2006.217 2006.226 2006.235
#> [1233] 2006.244 2006.253 2006.262 2006.271 2006.281 2006.290 2006.299 2006.308
#> [1241] 2006.317 2006.326 2006.335 2006.344 2006.353 2006.362 2006.371 2006.380
#> [1249] 2006.389 2006.399 2006.408 2006.417 2006.426 2006.436 2006.445 2006.454
#> [1257] 2006.463 2006.472 2006.482 2006.491 2006.500 2006.509 2006.518 2006.527
#> [1265] 2006.536 2006.545 2006.554 2006.563 2006.572 2006.581 2006.590 2006.599
#> [1273] 2006.608 2006.617 2006.626 2006.635 2006.644 2006.653 2006.661 2006.669
#> [1281] 2006.677 2006.684 2006.694 2006.703 2006.713 2006.722 2006.732 2006.741
#> [1289] 2006.751 2006.760 2006.769 2006.779 2006.788 2006.798 2006.807 2006.817
#> [1297] 2006.826 2006.835 2006.844 2006.853 2006.862 2006.871 2006.880 2006.889
#> [1305] 2006.898 2006.906 2006.915 2006.924 2006.933 2006.942 2006.951 2006.960
#> [1313] 2006.968 2006.976 2006.985 2006.994 2007.004 2007.013 2007.022 2007.032
#> [1321] 2007.041 2007.050 2007.060 2007.069 2007.079 2007.088 2007.097 2007.107
#> [1329] 2007.116 2007.126 2007.135 2007.145 2007.154 2007.164 2007.173 2007.182
#> [1337] 2007.192 2007.201 2007.211 2007.220 2007.230 2007.239 2007.248 2007.258
#> [1345] 2007.266 2007.275 2007.283 2007.292 2007.301 2007.310 2007.319 2007.328
#> [1353] 2007.337 2007.346 2007.355 2007.364 2007.373 2007.381 2007.389 2007.398
#> [1361] 2007.406 2007.415 2007.423 2007.431 2007.440 2007.448 2007.457 2007.466
#> [1369] 2007.474 2007.483 2007.491 2007.500 2007.509 2007.518 2007.528 2007.537
#> [1377] 2007.546 2007.555 2007.564 2007.574 2007.583 2007.592 2007.601 2007.610
#> [1385] 2007.620 2007.629 2007.638 2007.647 2007.656 2007.666 2007.675 2007.684
#> [1393] 2007.693 2007.702 2007.712 2007.721 2007.730 2007.739 2007.748 2007.758
#> [1401] 2007.767 2007.776 2007.785 2007.794 2007.804 2007.813 2007.822 2007.831
#> [1409] 2007.840 2007.849 2007.858 2007.867 2007.876 2007.885 2007.894 2007.903
#> [1417] 2007.912 2007.921 2007.930 2007.938 2007.947 2007.955 2007.965 2007.974
#> [1425] 2007.983 2007.993 2008.002 2008.011 2008.021 2008.029 2008.038 2008.047
#> [1433] 2008.056 2008.064 2008.073 2008.082 2008.091 2008.099 2008.108 2008.117
#> [1441] 2008.126 2008.135 2008.144 2008.153 2008.162 2008.171 2008.180 2008.189
#> [1449] 2008.198 2008.207 2008.216 2008.225 2008.234 2008.243 2008.252 2008.261
#> [1457] 2008.270 2008.279 2008.288 2008.297 2008.306 2008.314 2008.323 2008.333
#> [1465] 2008.342 2008.351 2008.360 2008.369 2008.378 2008.387 2008.397 2008.406
#> [1473] 2008.415 2008.424 2008.433 2008.442 2008.451 2008.460 2008.468 2008.476
#> [1481] 2008.485 2008.492 2008.500 2008.509 2008.519 2008.528 2008.537 2008.547
#> [1489] 2008.556 2008.565 2008.573 2008.582 2008.591 2008.599 2008.608 2008.617
#> [1497] 2008.625 2008.635 2008.644 2008.653 2008.662 2008.671 2008.680 2008.689
#> [1505] 2008.699 2008.708 2008.717 2008.726 2008.736 2008.745 2008.755 2008.764
#> [1513] 2008.774 2008.783 2008.792 2008.802 2008.811 2008.821 2008.830 2008.838
#> [1521] 2008.847 2008.856 2008.864 2008.873 2008.882 2008.890 2008.899 2008.908
#> [1529] 2008.917 2008.925 2008.933 2008.941 2008.949 2008.958 2008.966 2008.974
#> [1537] 2008.983 2008.992 2009.001 2009.010 2009.018 2009.027 2009.035 2009.044
#> [1545] 2009.052 2009.061 2009.070 2009.079 2009.087 2009.096 2009.105 2009.114
#> [1553] 2009.123 2009.132 2009.141 2009.150 2009.159 2009.168 2009.176 2009.185
#> [1561] 2009.193 2009.201 2009.209 2009.217 2009.226 2009.236 2009.245 2009.255
#> [1569] 2009.264 2009.272 2009.281 2009.290 2009.299 2009.308 2009.316 2009.325
#> [1577] 2009.333 2009.341 2009.349 2009.355 2009.364 2009.373 2009.381 2009.390
#> [1585] 2009.399 2009.407 2009.416 2009.425 2009.433 2009.442 2009.450 2009.458
#> [1593] 2009.467 2009.475 2009.483 2009.492 2009.500 2009.505 2009.514 2009.523
#> [1601] 2009.533 2009.542 2009.551 2009.561 2009.570 2009.579 2009.588 2009.598
#> [1609] 2009.607 2009.616 2009.626 2009.635 2009.644 2009.653 2009.663 2009.672
#> [1617] 2009.681 2009.691 2009.700 2009.709 2009.717 2009.726 2009.735 2009.744
#> [1625] 2009.752 2009.761 2009.770 2009.779 2009.788 2009.796 2009.806 2009.815
#> [1633] 2009.825 2009.834 2009.844 2009.853 2009.863 2009.872 2009.882 2009.889
#> [1641] 2009.897 2009.905 2009.909 2009.918 2009.928 2009.937 2009.946 2009.947
#> [1649] 2009.954 2009.962 2009.970 2009.977 2009.985 2009.993 2010.001 2010.008
#> [1657] 2010.016 2010.026 2010.035 2010.045 2010.054 2010.064 2010.073 2010.083
#> [1665] 2010.092 2010.102 2010.111 2010.120 2010.130 2010.138 2010.146 2010.153
#> [1673] 2010.162 2010.170 2010.179 2010.187 2010.195 2010.204 2010.212 2010.221
#> [1681] 2010.230 2010.240 2010.249 2010.258 2010.267 2010.276 2010.285 2010.295
#> [1689] 2010.304 2010.313 2010.322 2010.331 2010.340 2010.350 2010.359 2010.368
#> [1697] 2010.377 2010.387 2010.396 2010.405 2010.415 2010.424 2010.434 2010.443
#> [1705] 2010.453 2010.462 2010.472 2010.481 2010.491 2010.500 2010.509 2010.519
#> [1713] 2010.528 2010.538 2010.547 2010.557 2010.566 2010.573 2010.580 2010.588
#> [1721] 2010.594 2010.600 2010.609 2010.619 2010.628 2010.637 2010.647 2010.656
#> [1729] 2010.665 2010.674 2010.684 2010.693 2010.702 2010.712 2010.721 2010.730
#> [1737] 2010.737 2010.743 2010.752 2010.761 2010.770 2010.779 2010.788 2010.797
#> [1745] 2010.806 2010.815 2010.824 2010.833 2010.842 2010.851 2010.860 2010.869
#> [1753] 2010.877 2010.886 2010.895 2010.904 2010.913 2010.921 2010.929 2010.937
#> [1761] 2010.944 2010.952 2010.961 2010.969 2010.977 2010.986 2010.994 2011.002
#> [1769] 2011.011 2011.020 2011.028 2011.037 2011.046 2011.055 2011.063 2011.072
#> [1777] 2011.081 2011.090 2011.099 2011.107 2011.116 2011.125 2011.134 2011.143
#> [1785] 2011.152 2011.161 2011.170 2011.180 2011.189 2011.198 2011.208 2011.217
#> [1793] 2011.227 2011.236 2011.245 2011.255 2011.264 2011.273 2011.283 2011.292
#> [1801] 2011.300 2011.308 2011.316 2011.324 2011.332 2011.340 2011.348 2011.353
#> [1809] 2011.359 2011.369 2011.378 2011.387 2011.397 2011.406 2011.415 2011.425
#> [1817] 2011.434 2011.444 2011.453 2011.462 2011.472 2011.481 2011.491 2011.500
#> [1825] 2011.509 2011.519 2011.528 2011.537 2011.547 2011.556 2011.566 2011.575
#> [1833] 2011.584 2011.594 2011.603 2011.612 2011.620 2011.628 2011.636 2011.643
#> [1841] 2011.652 2011.660 2011.668 2011.675 2011.684 2011.694 2011.703 2011.713
#> [1849] 2011.722 2011.732 2011.741 2011.751 2011.760 2011.769 2011.779 2011.788
#> [1857] 2011.798 2011.807 2011.817 2011.826 2011.836 2011.845 2011.855 2011.864
#> [1865] 2011.874 2011.883 2011.892 2011.902 2011.911 2011.920 2011.928 2011.936
#> [1873] 2011.945 2011.954 2011.963 2011.972 2011.981 2011.990 2011.999 2012.008
#> [1881] 2012.017 2012.026 2012.035 2012.044 2012.053 2012.062 2012.071 2012.080
#> [1889] 2012.089 2012.098 2012.107 2012.116 2012.125 2012.134 2012.143 2012.146
#> [1897] 2012.155 2012.164 2012.173 2012.182 2012.191 2012.200 2012.208 2012.217
#> [1905] 2012.226 2012.235 2012.244 2012.253 2012.262 2012.271 2012.280 2012.289
#> [1913] 2012.298 2012.308 2012.317 2012.326 2012.336 2012.345 2012.354 2012.364
#> [1921] 2012.373 2012.382 2012.392 2012.401 2012.410 2012.420 2012.429 2012.438
#> [1929] 2012.448 2012.457 2012.466 2012.474 2012.483 2012.491 2012.500 2012.508
#> [1937] 2012.516 2012.524 2012.532 2012.540 2012.550 2012.559 2012.568 2012.577
#> [1945] 2012.586 2012.596 2012.604 2012.613 2012.622 2012.631 2012.639 2012.648
#> [1953] 2012.657 2012.666 2012.675 2012.683 2012.692 2012.701 2012.710 2012.719
#> [1961] 2012.727 2012.736 2012.745 2012.754 2012.762 2012.771 2012.780 2012.789
#> [1969] 2012.798 2012.806 2012.815 2012.824 2012.833 2012.841 2012.850 2012.859
#> [1977] 2012.867 2012.876 2012.884 2012.893 2012.901 2012.910 2012.918 2012.927
#> [1985] 2012.936 2012.945 2012.955 2012.964 2012.973 2012.982 2012.992 2013.001
#> [1993] 2013.010 2013.019 2013.029 2013.038 2013.047 2013.056 2013.066 2013.075
#> [2001] 2013.084 2013.093 2013.103 2013.112 2013.121 2013.130 2013.140 2013.149
#> [2009] 2013.158 2013.167 2013.177 2013.186 2013.195 2013.202 2013.211 2013.220
#> [2017] 2013.229 2013.238 2013.247 2013.256 2013.265 2013.274 2013.283 2013.292
#> [2025] 2013.301 2013.310 2013.319 2013.328 2013.337 2013.346 2013.355 2013.364
#> [2033] 2013.373 2013.382 2013.391 2013.400 2013.409 2013.418 2013.427 2013.436
#> [2041] 2013.446 2013.455 2013.464 2013.473 2013.482 2013.491
s5b <- plot_ppcheck_At(out, 2)
#>  [1]  0.5  1.5  2.5  3.5  4.5  5.5  6.5  7.5  8.5  9.5 10.5 11.5
#>  [1]  1  3 26 16 25  4 17  6 17 29 12  8
#> [1] 2014
#>    [1] 1995.000 1995.009 1995.019 1995.028 1995.038 1995.047 1995.056 1995.066
#>    [9] 1995.075 1995.085 1995.094 1995.103 1995.113 1995.122 1995.131 1995.141
#>   [17] 1995.150 1995.160 1995.169 1995.178 1995.188 1995.197 1995.207 1995.216
#>   [25] 1995.225 1995.235 1995.244 1995.254 1995.263 1995.272 1995.282 1995.291
#>   [33] 1995.300 1995.310 1995.319 1995.329 1995.338 1995.347 1995.357 1995.366
#>   [41] 1995.376 1995.385 1995.394 1995.404 1995.413 1995.423 1995.432 1995.441
#>   [49] 1995.451 1995.460 1995.469 1995.479 1995.488 1995.498 1995.507 1995.516
#>   [57] 1995.526 1995.535 1995.544 1995.553 1995.562 1995.572 1995.581 1995.590
#>   [65] 1995.599 1995.608 1995.617 1995.626 1995.635 1995.644 1995.653 1995.663
#>   [73] 1995.672 1995.681 1995.690 1995.699 1995.708 1995.717 1995.727 1995.736
#>   [81] 1995.745 1995.755 1995.764 1995.773 1995.783 1995.792 1995.801 1995.811
#>   [89] 1995.820 1995.829 1995.838 1995.848 1995.857 1995.866 1995.876 1995.885
#>   [97] 1995.894 1995.904 1995.913 1995.922 1995.932 1995.941 1995.950 1995.960
#>  [105] 1995.969 1995.978 1995.988 1995.997 1996.006 1996.015 1996.025 1996.034
#>  [113] 1996.044 1996.053 1996.063 1996.072 1996.081 1996.091 1996.100 1996.110
#>  [121] 1996.119 1996.129 1996.138 1996.147 1996.157 1996.166 1996.176 1996.185
#>  [129] 1996.195 1996.204 1996.213 1996.223 1996.232 1996.242 1996.251 1996.261
#>  [137] 1996.270 1996.279 1996.289 1996.298 1996.308 1996.317 1996.327 1996.336
#>  [145] 1996.345 1996.355 1996.364 1996.374 1996.383 1996.393 1996.402 1996.411
#>  [153] 1996.421 1996.430 1996.440 1996.449 1996.459 1996.468 1996.477 1996.487
#>  [161] 1996.496 1996.506 1996.515 1996.525 1996.534 1996.543 1996.553 1996.562
#>  [169] 1996.572 1996.581 1996.591 1996.600 1996.609 1996.619 1996.628 1996.638
#>  [177] 1996.647 1996.657 1996.666 1996.675 1996.685 1996.694 1996.704 1996.713
#>  [185] 1996.722 1996.731 1996.740 1996.749 1996.758 1996.768 1996.777 1996.786
#>  [193] 1996.796 1996.805 1996.815 1996.824 1996.834 1996.843 1996.853 1996.862
#>  [201] 1996.872 1996.881 1996.891 1996.900 1996.910 1996.919 1996.929 1996.938
#>  [209] 1996.948 1996.957 1996.967 1996.976 1996.986 1996.995 1997.005 1997.014
#>  [217] 1997.024 1997.033 1997.043 1997.052 1997.062 1997.071 1997.081 1997.090
#>  [225] 1997.100 1997.109 1997.119 1997.128 1997.137 1997.147 1997.156 1997.165
#>  [233] 1997.175 1997.184 1997.193 1997.203 1997.212 1997.221 1997.231 1997.240
#>  [241] 1997.249 1997.259 1997.268 1997.277 1997.287 1997.296 1997.305 1997.315
#>  [249] 1997.324 1997.334 1997.343 1997.352 1997.362 1997.371 1997.380 1997.390
#>  [257] 1997.399 1997.408 1997.418 1997.427 1997.436 1997.446 1997.455 1997.464
#>  [265] 1997.474 1997.483 1997.492 1997.501 1997.510 1997.519 1997.528 1997.537
#>  [273] 1997.546 1997.555 1997.564 1997.573 1997.582 1997.591 1997.600 1997.609
#>  [281] 1997.618 1997.627 1997.636 1997.645 1997.654 1997.664 1997.673 1997.682
#>  [289] 1997.692 1997.701 1997.710 1997.719 1997.729 1997.738 1997.747 1997.757
#>  [297] 1997.766 1997.775 1997.785 1997.794 1997.803 1997.813 1997.822 1997.831
#>  [305] 1997.840 1997.850 1997.859 1997.868 1997.878 1997.887 1997.896 1997.906
#>  [313] 1997.915 1997.924 1997.934 1997.943 1997.952 1997.962 1997.971 1997.981
#>  [321] 1997.990 1998.000 1998.009 1998.019 1998.028 1998.038 1998.047 1998.057
#>  [329] 1998.066 1998.076 1998.085 1998.095 1998.104 1998.114 1998.123 1998.133
#>  [337] 1998.142 1998.151 1998.161 1998.170 1998.180 1998.189 1998.199 1998.208
#>  [345] 1998.218 1998.227 1998.237 1998.246 1998.256 1998.265 1998.275 1998.284
#>  [353] 1998.294 1998.303 1998.313 1998.322 1998.332 1998.341 1998.351 1998.360
#>  [361] 1998.370 1998.379 1998.389 1998.398 1998.408 1998.417 1998.427 1998.436
#>  [369] 1998.446 1998.455 1998.464 1998.474 1998.483 1998.493 1998.502 1998.512
#>  [377] 1998.521 1998.531 1998.540 1998.550 1998.559 1998.569 1998.578 1998.588
#>  [385] 1998.597 1998.607 1998.616 1998.626 1998.635 1998.639 1998.648 1998.657
#>  [393] 1998.667 1998.676 1998.685 1998.694 1998.704 1998.713 1998.722 1998.732
#>  [401] 1998.741 1998.751 1998.760 1998.769 1998.778 1998.787 1998.795 1998.804
#>  [409] 1998.813 1998.822 1998.831 1998.840 1998.849 1998.859 1998.868 1998.877
#>  [417] 1998.887 1998.896 1998.905 1998.915 1998.924 1998.934 1998.943 1998.952
#>  [425] 1998.962 1998.971 1998.980 1998.990 1998.999 1999.008 1999.018 1999.027
#>  [433] 1999.037 1999.046 1999.055 1999.065 1999.074 1999.083 1999.093 1999.102
#>  [441] 1999.112 1999.121 1999.130 1999.140 1999.149 1999.158 1999.168 1999.177
#>  [449] 1999.186 1999.196 1999.205 1999.215 1999.224 1999.233 1999.243 1999.252
#>  [457] 1999.261 1999.271 1999.280 1999.289 1999.299 1999.307 1999.314 1999.322
#>  [465] 1999.330 1999.338 1999.347 1999.356 1999.365 1999.375 1999.384 1999.393
#>  [473] 1999.402 1999.411 1999.420 1999.429 1999.439 1999.448 1999.457 1999.466
#>  [481] 1999.475 1999.484 1999.493 1999.503 1999.512 1999.521 1999.530 1999.538
#>  [489] 1999.547 1999.556 1999.565 1999.574 1999.582 1999.589 1999.596 1999.605
#>  [497] 1999.614 1999.624 1999.633 1999.642 1999.651 1999.661 1999.670 1999.679
#>  [505] 1999.689 1999.698 1999.707 1999.716 1999.726 1999.735 1999.744 1999.754
#>  [513] 1999.763 1999.772 1999.781 1999.791 1999.800 1999.809 1999.819 1999.828
#>  [521] 1999.837 1999.846 1999.855 1999.864 1999.873 1999.882 1999.891 1999.900
#>  [529] 1999.909 1999.918 1999.927 1999.936 1999.945 1999.954 1999.963 1999.972
#>  [537] 1999.981 1999.988 1999.995 2000.004 2000.014 2000.023 2000.032 2000.042
#>  [545] 2000.051 2000.061 2000.070 2000.079 2000.089 2000.098 2000.108 2000.117
#>  [553] 2000.126 2000.136 2000.145 2000.155 2000.164 2000.173 2000.183 2000.192
#>  [561] 2000.201 2000.211 2000.220 2000.230 2000.238 2000.247 2000.256 2000.265
#>  [569] 2000.273 2000.281 2000.289 2000.297 2000.304 2000.312 2000.322 2000.331
#>  [577] 2000.341 2000.350 2000.359 2000.369 2000.378 2000.388 2000.397 2000.406
#>  [585] 2000.416 2000.425 2000.434 2000.444 2000.453 2000.463 2000.472 2000.481
#>  [593] 2000.490 2000.500 2000.509 2000.518 2000.527 2000.536 2000.545 2000.554
#>  [601] 2000.563 2000.573 2000.582 2000.591 2000.600 2000.609 2000.618 2000.627
#>  [609] 2000.636 2000.646 2000.655 2000.664 2000.673 2000.682 2000.691 2000.700
#>  [617] 2000.709 2000.718 2000.724 2000.731 2000.737 2000.746 2000.755 2000.763
#>  [625] 2000.772 2000.781 2000.790 2000.798 2000.807 2000.816 2000.825 2000.835
#>  [633] 2000.844 2000.854 2000.863 2000.872 2000.882 2000.891 2000.901 2000.910
#>  [641] 2000.920 2000.929 2000.939 2000.948 2000.958 2000.967 2000.977 2000.986
#>  [649] 2000.995 2001.005 2001.014 2001.023 2001.033 2001.042 2001.051 2001.061
#>  [657] 2001.070 2001.079 2001.089 2001.098 2001.108 2001.117 2001.126 2001.135
#>  [665] 2001.144 2001.154 2001.163 2001.172 2001.181 2001.190 2001.199 2001.208
#>  [673] 2001.216 2001.225 2001.234 2001.243 2001.252 2001.261 2001.270 2001.279
#>  [681] 2001.288 2001.297 2001.306 2001.310 2001.319 2001.328 2001.337 2001.346
#>  [689] 2001.354 2001.363 2001.372 2001.381 2001.390 2001.399 2001.408 2001.418
#>  [697] 2001.427 2001.436 2001.446 2001.455 2001.464 2001.474 2001.483 2001.489
#>  [705] 2001.494 2001.502 2001.510 2001.518 2001.526 2001.535 2001.543 2001.549
#>  [713] 2001.555 2001.564 2001.573 2001.583 2001.592 2001.602 2001.611 2001.620
#>  [721] 2001.630 2001.639 2001.648 2001.658 2001.667 2001.677 2001.686 2001.695
#>  [729] 2001.705 2001.714 2001.723 2001.733 2001.741 2001.750 2001.760 2001.769
#>  [737] 2001.778 2001.787 2001.796 2001.806 2001.815 2001.824 2001.833 2001.842
#>  [745] 2001.851 2001.861 2001.870 2001.879 2001.888 2001.897 2001.907 2001.916
#>  [753] 2001.925 2001.932 2001.940 2001.948 2001.957 2001.966 2001.975 2001.984
#>  [761] 2001.993 2002.002 2002.011 2002.020 2002.030 2002.039 2002.047 2002.056
#>  [769] 2002.065 2002.074 2002.083 2002.092 2002.101 2002.110 2002.119 2002.128
#>  [777] 2002.138 2002.147 2002.157 2002.166 2002.175 2002.185 2002.194 2002.204
#>  [785] 2002.213 2002.223 2002.232 2002.242 2002.251 2002.261 2002.270 2002.280
#>  [793] 2002.289 2002.298 2002.308 2002.317 2002.327 2002.336 2002.346 2002.355
#>  [801] 2002.365 2002.367 2002.372 2002.381 2002.390 2002.399 2002.409 2002.418
#>  [809] 2002.427 2002.436 2002.445 2002.455 2002.464 2002.473 2002.482 2002.491
#>  [817] 2002.500 2002.505 2002.514 2002.522 2002.531 2002.540 2002.548 2002.557
#>  [825] 2002.565 2002.574 2002.583 2002.590 2002.597 2002.604 2002.613 2002.621
#>  [833] 2002.629 2002.638 2002.646 2002.654 2002.663 2002.671 2002.679 2002.688
#>  [841] 2002.696 2002.705 2002.712 2002.719 2002.726 2002.728 2002.737 2002.745
#>  [849] 2002.754 2002.762 2002.771 2002.781 2002.790 2002.799 2002.808 2002.817
#>  [857] 2002.827 2002.836 2002.845 2002.854 2002.863 2002.873 2002.882 2002.891
#>  [865] 2002.900 2002.907 2002.913 2002.917 2002.924 2002.932 2002.940 2002.947
#>  [873] 2002.956 2002.965 2002.973 2002.982 2002.991 2003.001 2003.010 2003.019
#>  [881] 2003.028 2003.037 2003.046 2003.056 2003.065 2003.074 2003.083 2003.092
#>  [889] 2003.102 2003.111 2003.120 2003.129 2003.138 2003.147 2003.157 2003.166
#>  [897] 2003.175 2003.184 2003.193 2003.200 2003.210 2003.219 2003.228 2003.238
#>  [905] 2003.247 2003.256 2003.265 2003.275 2003.284 2003.293 2003.302 2003.312
#>  [913] 2003.321 2003.330 2003.339 2003.349 2003.358 2003.367 2003.376 2003.385
#>  [921] 2003.395 2003.404 2003.413 2003.422 2003.430 2003.439 2003.448 2003.457
#>  [929] 2003.466 2003.475 2003.484 2003.492 2003.500 2003.508 2003.517 2003.525
#>  [937] 2003.530 2003.535 2003.544 2003.553 2003.558 2003.567 2003.576 2003.585
#>  [945] 2003.594 2003.603 2003.612 2003.621 2003.631 2003.640 2003.649 2003.658
#>  [953] 2003.667 2003.676 2003.685 2003.693 2003.701 2003.708 2003.716 2003.717
#>  [961] 2003.726 2003.735 2003.744 2003.753 2003.762 2003.771 2003.780 2003.789
#>  [969] 2003.798 2003.807 2003.817 2003.826 2003.835 2003.844 2003.854 2003.863
#>  [977] 2003.872 2003.881 2003.891 2003.900 2003.909 2003.918 2003.928 2003.937
#>  [985] 2003.946 2003.951 2003.960 2003.969 2003.978 2003.987 2003.993 2004.001
#>  [993] 2004.010 2004.018 2004.026 2004.034 2004.043 2004.052 2004.061 2004.071
#> [1001] 2004.080 2004.090 2004.099 2004.108 2004.118 2004.127 2004.136 2004.146
#> [1009] 2004.155 2004.164 2004.174 2004.183 2004.192 2004.202 2004.211 2004.220
#> [1017] 2004.230 2004.239 2004.248 2004.257 2004.266 2004.275 2004.284 2004.293
#> [1025] 2004.301 2004.310 2004.319 2004.327 2004.337 2004.346 2004.356 2004.365
#> [1033] 2004.375 2004.384 2004.393 2004.402 2004.411 2004.419 2004.428 2004.437
#> [1041] 2004.445 2004.454 2004.463 2004.472 2004.480 2004.489 2004.498 2004.500
#> [1049] 2004.502 2004.511 2004.520 2004.529 2004.536 2004.543 2004.552 2004.560
#> [1057] 2004.569 2004.578 2004.587 2004.596 2004.605 2004.614 2004.623 2004.632
#> [1065] 2004.641 2004.650 2004.658 2004.667 2004.676 2004.686 2004.695 2004.704
#> [1073] 2004.713 2004.723 2004.732 2004.741 2004.751 2004.760 2004.769 2004.778
#> [1081] 2004.788 2004.797 2004.806 2004.816 2004.825 2004.834 2004.843 2004.853
#> [1089] 2004.862 2004.871 2004.881 2004.890 2004.899 2004.909 2004.918 2004.927
#> [1097] 2004.936 2004.946 2004.955 2004.965 2004.974 2004.984 2004.993 2005.003
#> [1105] 2005.012 2005.022 2005.031 2005.041 2005.050 2005.060 2005.069 2005.078
#> [1113] 2005.088 2005.097 2005.107 2005.116 2005.126 2005.135 2005.144 2005.153
#> [1121] 2005.162 2005.171 2005.180 2005.189 2005.198 2005.208 2005.217 2005.226
#> [1129] 2005.235 2005.244 2005.253 2005.262 2005.271 2005.280 2005.289 2005.298
#> [1137] 2005.307 2005.316 2005.325 2005.334 2005.344 2005.353 2005.362 2005.371
#> [1145] 2005.380 2005.390 2005.399 2005.408 2005.417 2005.426 2005.436 2005.445
#> [1153] 2005.454 2005.463 2005.472 2005.482 2005.491 2005.500 2005.509 2005.518
#> [1161] 2005.527 2005.536 2005.545 2005.554 2005.562 2005.571 2005.578 2005.585
#> [1169] 2005.592 2005.601 2005.611 2005.620 2005.630 2005.639 2005.649 2005.658
#> [1177] 2005.668 2005.677 2005.687 2005.696 2005.706 2005.715 2005.725 2005.734
#> [1185] 2005.744 2005.753 2005.763 2005.772 2005.782 2005.791 2005.801 2005.810
#> [1193] 2005.820 2005.829 2005.839 2005.848 2005.858 2005.867 2005.877 2005.885
#> [1201] 2005.894 2005.902 2005.911 2005.920 2005.929 2005.938 2005.947 2005.956
#> [1209] 2005.965 2005.974 2005.983 2005.992 2006.002 2006.011 2006.020 2006.028
#> [1217] 2006.036 2006.044 2006.049 2006.054 2006.063 2006.073 2006.082 2006.091
#> [1225] 2006.100 2006.109 2006.119 2006.128 2006.137 2006.146 2006.156 2006.165
#> [1233] 2006.174 2006.183 2006.193 2006.202 2006.211 2006.220 2006.229 2006.239
#> [1241] 2006.248 2006.257 2006.266 2006.276 2006.283 2006.291 2006.299 2006.308
#> [1249] 2006.317 2006.326 2006.334 2006.343 2006.352 2006.361 2006.370 2006.379
#> [1257] 2006.387 2006.396 2006.405 2006.413 2006.421 2006.428 2006.437 2006.446
#> [1265] 2006.455 2006.464 2006.473 2006.482 2006.491 2006.500 2006.509 2006.519
#> [1273] 2006.528 2006.538 2006.547 2006.556 2006.566 2006.575 2006.585 2006.594
#> [1281] 2006.603 2006.613 2006.622 2006.632 2006.637 2006.642 2006.651 2006.661
#> [1289] 2006.670 2006.679 2006.688 2006.698 2006.707 2006.716 2006.726 2006.735
#> [1297] 2006.744 2006.753 2006.763 2006.772 2006.781 2006.791 2006.800 2006.810
#> [1305] 2006.819 2006.828 2006.838 2006.847 2006.857 2006.866 2006.875 2006.885
#> [1313] 2006.894 2006.904 2006.913 2006.923 2006.932 2006.941 2006.951 2006.960
#> [1321] 2006.970 2006.979 2006.988 2006.998 2007.007 2007.017 2007.026 2007.036
#> [1329] 2007.045 2007.054 2007.064 2007.073 2007.083 2007.092 2007.101 2007.110
#> [1337] 2007.119 2007.127 2007.136 2007.144 2007.153 2007.161 2007.170 2007.179
#> [1345] 2007.187 2007.190 2007.199 2007.208 2007.217 2007.225 2007.234 2007.241
#> [1353] 2007.247 2007.256 2007.265 2007.274 2007.284 2007.293 2007.302 2007.311
#> [1361] 2007.320 2007.329 2007.338 2007.348 2007.357 2007.366 2007.375 2007.384
#> [1369] 2007.393 2007.402 2007.412 2007.421 2007.430 2007.438 2007.446 2007.454
#> [1377] 2007.462 2007.470 2007.477 2007.485 2007.492 2007.500 2007.509 2007.519
#> [1385] 2007.528 2007.538 2007.547 2007.557 2007.566 2007.575 2007.585 2007.594
#> [1393] 2007.604 2007.613 2007.623 2007.632 2007.641 2007.651 2007.660 2007.670
#> [1401] 2007.679 2007.689 2007.698 2007.707 2007.717 2007.726 2007.735 2007.744
#> [1409] 2007.753 2007.762 2007.771 2007.779 2007.788 2007.797 2007.806 2007.815
#> [1417] 2007.824 2007.833 2007.841 2007.849 2007.858 2007.866 2007.874 2007.882
#> [1425] 2007.891 2007.900 2007.910 2007.919 2007.928 2007.937 2007.946 2007.956
#> [1433] 2007.965 2007.974 2007.983 2007.992 2008.002 2008.011 2008.020 2008.029
#> [1441] 2008.039 2008.048 2008.058 2008.067 2008.076 2008.086 2008.095 2008.105
#> [1449] 2008.114 2008.123 2008.133 2008.142 2008.152 2008.161 2008.170 2008.180
#> [1457] 2008.189 2008.199 2008.207 2008.215 2008.224 2008.232 2008.240 2008.249
#> [1465] 2008.257 2008.265 2008.274 2008.283 2008.290 2008.299 2008.308 2008.317
#> [1473] 2008.326 2008.336 2008.345 2008.354 2008.363 2008.372 2008.381 2008.390
#> [1481] 2008.400 2008.409 2008.418 2008.427 2008.436 2008.445 2008.454 2008.463
#> [1489] 2008.473 2008.482 2008.491 2008.500 2008.505 2008.510 2008.519 2008.528
#> [1497] 2008.537 2008.546 2008.555 2008.564 2008.573 2008.582 2008.591 2008.600
#> [1505] 2008.609 2008.618 2008.626 2008.635 2008.643 2008.653 2008.662 2008.671
#> [1513] 2008.680 2008.689 2008.699 2008.708 2008.717 2008.726 2008.735 2008.744
#> [1521] 2008.752 2008.760 2008.767 2008.775 2008.783 2008.792 2008.802 2008.811
#> [1529] 2008.820 2008.830 2008.839 2008.849 2008.858 2008.868 2008.877 2008.887
#> [1537] 2008.896 2008.905 2008.915 2008.924 2008.933 2008.942 2008.950 2008.959
#> [1545] 2008.968 2008.977 2008.985 2008.994 2009.003 2009.012 2009.020 2009.029
#> [1553] 2009.038 2009.047 2009.055 2009.064 2009.072 2009.081 2009.089 2009.097
#> [1561] 2009.106 2009.114 2009.122 2009.132 2009.141 2009.150 2009.159 2009.168
#> [1569] 2009.178 2009.187 2009.196 2009.205 2009.215 2009.224 2009.233 2009.242
#> [1577] 2009.252 2009.261 2009.270 2009.279 2009.289 2009.298 2009.307 2009.316
#> [1585] 2009.326 2009.335 2009.344 2009.353 2009.361 2009.369 2009.377 2009.386
#> [1593] 2009.395 2009.404 2009.413 2009.422 2009.431 2009.440 2009.449 2009.458
#> [1601] 2009.467 2009.475 2009.484 2009.492 2009.500 2009.509 2009.519 2009.528
#> [1609] 2009.538 2009.547 2009.557 2009.566 2009.576 2009.585 2009.594 2009.604
#> [1617] 2009.613 2009.623 2009.632 2009.642 2009.651 2009.661 2009.670 2009.679
#> [1625] 2009.689 2009.698 2009.708 2009.717 2009.727 2009.736 2009.746 2009.755
#> [1633] 2009.764 2009.774 2009.783 2009.792 2009.802 2009.811 2009.820 2009.830
#> [1641] 2009.839 2009.848 2009.858 2009.867 2009.877 2009.886 2009.895 2009.905
#> [1649] 2009.914 2009.923 2009.933 2009.942 2009.951 2009.961 2009.970 2009.980
#> [1657] 2009.989 2009.998 2010.002 2010.008 2010.014 2010.021 2010.029 2010.035
#> [1665] 2010.042 2010.048 2010.058 2010.067 2010.077 2010.086 2010.096 2010.105
#> [1673] 2010.114 2010.124 2010.133 2010.143 2010.152 2010.162 2010.171 2010.181
#> [1681] 2010.190 2010.200 2010.209 2010.219 2010.228 2010.237 2010.247 2010.256
#> [1689] 2010.266 2010.275 2010.285 2010.294 2010.304 2010.313 2010.323 2010.332
#> [1697] 2010.342 2010.351 2010.360 2010.370 2010.374 2010.377 2010.386 2010.395
#> [1705] 2010.404 2010.412 2010.421 2010.430 2010.439 2010.447 2010.456 2010.465
#> [1713] 2010.474 2010.482 2010.491 2010.500 2010.508 2010.517 2010.525 2010.534
#> [1721] 2010.542 2010.552 2010.561 2010.570 2010.580 2010.589 2010.598 2010.608
#> [1729] 2010.617 2010.626 2010.636 2010.645 2010.654 2010.664 2010.673 2010.683
#> [1737] 2010.692 2010.701 2010.711 2010.720 2010.729 2010.739 2010.748 2010.757
#> [1745] 2010.767 2010.776 2010.785 2010.795 2010.804 2010.814 2010.823 2010.832
#> [1753] 2010.841 2010.850 2010.859 2010.868 2010.877 2010.887 2010.896 2010.905
#> [1761] 2010.914 2010.923 2010.932 2010.942 2010.951 2010.960 2010.969 2010.978
#> [1769] 2010.987 2010.997 2011.006 2011.015 2011.024 2011.033 2011.043 2011.052
#> [1777] 2011.061 2011.070 2011.079 2011.089 2011.097 2011.105 2011.113 2011.122
#> [1785] 2011.131 2011.141 2011.150 2011.160 2011.169 2011.179 2011.188 2011.198
#> [1793] 2011.207 2011.216 2011.226 2011.235 2011.245 2011.254 2011.264 2011.273
#> [1801] 2011.283 2011.292 2011.302 2011.311 2011.320 2011.330 2011.339 2011.349
#> [1809] 2011.358 2011.368 2011.377 2011.387 2011.396 2011.405 2011.415 2011.424
#> [1817] 2011.434 2011.443 2011.453 2011.462 2011.472 2011.481 2011.491 2011.500
#> [1825] 2011.509 2011.519 2011.528 2011.538 2011.547 2011.557 2011.566 2011.575
#> [1833] 2011.585 2011.594 2011.604 2011.613 2011.623 2011.632 2011.642 2011.651
#> [1841] 2011.660 2011.670 2011.679 2011.689 2011.698 2011.708 2011.717 2011.726
#> [1849] 2011.736 2011.745 2011.755 2011.764 2011.774 2011.783 2011.792 2011.802
#> [1857] 2011.811 2011.821 2011.830 2011.840 2011.849 2011.858 2011.868 2011.877
#> [1865] 2011.887 2011.896 2011.906 2011.915 2011.925 2011.934 2011.943 2011.953
#> [1873] 2011.962 2011.972 2011.981 2011.991 2012.000 2012.009 2012.019 2012.028
#> [1881] 2012.038 2012.047 2012.057 2012.066 2012.075 2012.085 2012.094 2012.104
#> [1889] 2012.113 2012.123 2012.132 2012.142 2012.151 2012.160 2012.170 2012.179
#> [1897] 2012.189 2012.198 2012.208 2012.217 2012.226 2012.236 2012.245 2012.255
#> [1905] 2012.264 2012.274 2012.283 2012.292 2012.302 2012.311 2012.321 2012.330
#> [1913] 2012.340 2012.349 2012.358 2012.368 2012.377 2012.387 2012.396 2012.406
#> [1921] 2012.415 2012.425 2012.434 2012.443 2012.453 2012.462 2012.472 2012.481
#> [1929] 2012.491 2012.500 2012.509 2012.519 2012.528 2012.538 2012.547 2012.557
#> [1937] 2012.566 2012.575 2012.585 2012.594 2012.604 2012.613 2012.623 2012.632
#> [1945] 2012.642 2012.651 2012.660 2012.670 2012.679 2012.689 2012.698 2012.708
#> [1953] 2012.717 2012.726 2012.736 2012.745 2012.755 2012.764 2012.774 2012.783
#> [1961] 2012.792 2012.802 2012.811 2012.821 2012.830 2012.840 2012.849 2012.858
#> [1969] 2012.868 2012.877 2012.887 2012.896 2012.906 2012.915 2012.925 2012.934
#> [1977] 2012.943 2012.953 2012.962 2012.972 2012.981 2012.991 2013.000 2013.009
#> [1985] 2013.019 2013.028 2013.038 2013.047 2013.057 2013.066 2013.075 2013.085
#> [1993] 2013.094 2013.104 2013.113 2013.123 2013.132 2013.142 2013.151 2013.160
#> [2001] 2013.170 2013.179 2013.189 2013.198 2013.208 2013.217 2013.226 2013.236
#> [2009] 2013.245 2013.255 2013.264 2013.274 2013.283 2013.292 2013.302 2013.311
#> [2017] 2013.321 2013.330 2013.340 2013.349 2013.358 2013.368 2013.377 2013.387
#> [2025] 2013.396 2013.406 2013.415 2013.425 2013.434 2013.443 2013.453 2013.462
#> [2033] 2013.472 2013.481 2013.491
s5c <- plot_ppcheck_At(out, 3)
#> [1] 0.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5
#> [1]  5  8  5 13  3 13  7  8  3
#> [1] 2014
#>    [1] 1995.000 1995.009 1995.019 1995.028 1995.038 1995.047 1995.057 1995.066
#>    [9] 1995.076 1995.085 1995.095 1995.104 1995.114 1995.123 1995.133 1995.142
#>   [17] 1995.152 1995.161 1995.171 1995.180 1995.190 1995.199 1995.209 1995.218
#>   [25] 1995.228 1995.237 1995.247 1995.256 1995.266 1995.275 1995.285 1995.294
#>   [33] 1995.304 1995.313 1995.323 1995.332 1995.342 1995.351 1995.361 1995.370
#>   [41] 1995.380 1995.389 1995.399 1995.408 1995.418 1995.427 1995.437 1995.446
#>   [49] 1995.456 1995.465 1995.475 1995.484 1995.494 1995.503 1995.513 1995.522
#>   [57] 1995.532 1995.541 1995.551 1995.560 1995.570 1995.579 1995.589 1995.598
#>   [65] 1995.608 1995.617 1995.627 1995.636 1995.646 1995.655 1995.665 1995.674
#>   [73] 1995.684 1995.693 1995.703 1995.712 1995.722 1995.731 1995.741 1995.750
#>   [81] 1995.760 1995.769 1995.779 1995.788 1995.798 1995.807 1995.817 1995.826
#>   [89] 1995.836 1995.845 1995.855 1995.864 1995.874 1995.883 1995.893 1995.902
#>   [97] 1995.912 1995.921 1995.931 1995.940 1995.950 1995.959 1995.969 1995.978
#>  [105] 1995.988 1995.997 1996.007 1996.016 1996.026 1996.035 1996.045 1996.054
#>  [113] 1996.064 1996.073 1996.083 1996.092 1996.102 1996.111 1996.121 1996.130
#>  [121] 1996.140 1996.149 1996.159 1996.168 1996.178 1996.187 1996.197 1996.206
#>  [129] 1996.216 1996.225 1996.235 1996.244 1996.254 1996.263 1996.273 1996.282
#>  [137] 1996.292 1996.301 1996.311 1996.320 1996.330 1996.339 1996.349 1996.358
#>  [145] 1996.368 1996.377 1996.387 1996.396 1996.406 1996.415 1996.425 1996.434
#>  [153] 1996.444 1996.453 1996.463 1996.472 1996.482 1996.491 1996.501 1996.510
#>  [161] 1996.520 1996.529 1996.539 1996.548 1996.558 1996.567 1996.577 1996.586
#>  [169] 1996.596 1996.605 1996.615 1996.624 1996.634 1996.643 1996.653 1996.662
#>  [177] 1996.672 1996.681 1996.691 1996.700 1996.710 1996.719 1996.729 1996.738
#>  [185] 1996.748 1996.757 1996.767 1996.776 1996.786 1996.795 1996.805 1996.814
#>  [193] 1996.824 1996.833 1996.843 1996.852 1996.862 1996.871 1996.881 1996.890
#>  [201] 1996.900 1996.909 1996.919 1996.928 1996.938 1996.947 1996.957 1996.966
#>  [209] 1996.976 1996.985 1996.995 1997.004 1997.014 1997.023 1997.033 1997.042
#>  [217] 1997.052 1997.061 1997.071 1997.080 1997.090 1997.099 1997.109 1997.118
#>  [225] 1997.128 1997.137 1997.147 1997.156 1997.166 1997.175 1997.185 1997.194
#>  [233] 1997.204 1997.213 1997.223 1997.232 1997.242 1997.251 1997.261 1997.270
#>  [241] 1997.280 1997.289 1997.299 1997.308 1997.318 1997.327 1997.337 1997.346
#>  [249] 1997.356 1997.365 1997.375 1997.384 1997.394 1997.403 1997.413 1997.422
#>  [257] 1997.432 1997.441 1997.451 1997.460 1997.470 1997.479 1997.488 1997.498
#>  [265] 1997.507 1997.517 1997.526 1997.536 1997.545 1997.555 1997.564 1997.574
#>  [273] 1997.583 1997.593 1997.602 1997.612 1997.621 1997.631 1997.640 1997.650
#>  [281] 1997.659 1997.669 1997.678 1997.688 1997.697 1997.707 1997.716 1997.726
#>  [289] 1997.735 1997.745 1997.754 1997.764 1997.773 1997.783 1997.792 1997.801
#>  [297] 1997.811 1997.820 1997.830 1997.839 1997.849 1997.858 1997.868 1997.877
#>  [305] 1997.887 1997.896 1997.906 1997.915 1997.925 1997.934 1997.944 1997.953
#>  [313] 1997.963 1997.972 1997.982 1997.991 1998.001 1998.010 1998.019 1998.029
#>  [321] 1998.038 1998.048 1998.057 1998.067 1998.076 1998.086 1998.095 1998.105
#>  [329] 1998.114 1998.124 1998.133 1998.143 1998.152 1998.162 1998.171 1998.181
#>  [337] 1998.190 1998.200 1998.209 1998.219 1998.228 1998.237 1998.247 1998.256
#>  [345] 1998.266 1998.275 1998.285 1998.294 1998.304 1998.313 1998.323 1998.332
#>  [353] 1998.342 1998.351 1998.361 1998.370 1998.380 1998.389 1998.399 1998.408
#>  [361] 1998.418 1998.427 1998.437 1998.446 1998.455 1998.465 1998.474 1998.484
#>  [369] 1998.493 1998.503 1998.512 1998.522 1998.531 1998.541 1998.550 1998.560
#>  [377] 1998.569 1998.579 1998.588 1998.598 1998.607 1998.617 1998.626 1998.636
#>  [385] 1998.645 1998.655 1998.664 1998.674 1998.683 1998.692 1998.702 1998.711
#>  [393] 1998.721 1998.730 1998.740 1998.749 1998.759 1998.768 1998.778 1998.787
#>  [401] 1998.797 1998.806 1998.816 1998.825 1998.835 1998.844 1998.854 1998.863
#>  [409] 1998.873 1998.882 1998.892 1998.901 1998.910 1998.920 1998.929 1998.939
#>  [417] 1998.948 1998.958 1998.967 1998.977 1998.986 1998.996 1999.005 1999.015
#>  [425] 1999.024 1999.034 1999.043 1999.053 1999.062 1999.072 1999.081 1999.091
#>  [433] 1999.100 1999.110 1999.119 1999.128 1999.138 1999.147 1999.157 1999.166
#>  [441] 1999.176 1999.185 1999.195 1999.204 1999.214 1999.223 1999.233 1999.242
#>  [449] 1999.252 1999.261 1999.271 1999.280 1999.290 1999.299 1999.309 1999.318
#>  [457] 1999.328 1999.337 1999.346 1999.356 1999.365 1999.375 1999.384 1999.394
#>  [465] 1999.403 1999.413 1999.422 1999.432 1999.441 1999.451 1999.460 1999.470
#>  [473] 1999.479 1999.489 1999.498 1999.508 1999.517 1999.526 1999.536 1999.545
#>  [481] 1999.555 1999.564 1999.574 1999.583 1999.593 1999.602 1999.612 1999.621
#>  [489] 1999.631 1999.640 1999.650 1999.659 1999.669 1999.678 1999.687 1999.697
#>  [497] 1999.706 1999.716 1999.725 1999.735 1999.744 1999.754 1999.763 1999.773
#>  [505] 1999.782 1999.791 1999.801 1999.810 1999.819 1999.829 1999.838 1999.847
#>  [513] 1999.857 1999.866 1999.875 1999.885 1999.894 1999.904 1999.913 1999.922
#>  [521] 1999.932 1999.941 1999.950 1999.960 1999.969 1999.978 1999.988 1999.997
#>  [529] 2000.006 2000.016 2000.025 2000.034 2000.044 2000.053 2000.062 2000.072
#>  [537] 2000.081 2000.090 2000.100 2000.109 2000.118 2000.128 2000.137 2000.146
#>  [545] 2000.156 2000.165 2000.174 2000.184 2000.193 2000.203 2000.212 2000.221
#>  [553] 2000.231 2000.240 2000.249 2000.259 2000.268 2000.277 2000.287 2000.296
#>  [561] 2000.305 2000.315 2000.324 2000.333 2000.343 2000.352 2000.362 2000.371
#>  [569] 2000.380 2000.390 2000.399 2000.409 2000.418 2000.428 2000.437 2000.447
#>  [577] 2000.456 2000.466 2000.475 2000.484 2000.494 2000.503 2000.513 2000.522
#>  [585] 2000.532 2000.541 2000.551 2000.560 2000.569 2000.579 2000.588 2000.598
#>  [593] 2000.607 2000.617 2000.626 2000.636 2000.645 2000.655 2000.664 2000.673
#>  [601] 2000.683 2000.692 2000.702 2000.711 2000.721 2000.730 2000.740 2000.749
#>  [609] 2000.758 2000.768 2000.777 2000.787 2000.796 2000.806 2000.815 2000.825
#>  [617] 2000.834 2000.844 2000.853 2000.862 2000.872 2000.881 2000.891 2000.900
#>  [625] 2000.910 2000.919 2000.929 2000.938 2000.948 2000.957 2000.967 2000.976
#>  [633] 2000.985 2000.995 2001.004 2001.014 2001.023 2001.033 2001.042 2001.052
#>  [641] 2001.061 2001.071 2001.080 2001.090 2001.099 2001.109 2001.118 2001.128
#>  [649] 2001.137 2001.146 2001.156 2001.165 2001.175 2001.184 2001.194 2001.203
#>  [657] 2001.213 2001.222 2001.232 2001.241 2001.251 2001.260 2001.270 2001.279
#>  [665] 2001.289 2001.298 2001.307 2001.317 2001.326 2001.336 2001.345 2001.355
#>  [673] 2001.364 2001.374 2001.383 2001.393 2001.402 2001.412 2001.421 2001.431
#>  [681] 2001.440 2001.449 2001.459 2001.468 2001.478 2001.487 2001.497 2001.506
#>  [689] 2001.516 2001.525 2001.535 2001.544 2001.553 2001.563 2001.572 2001.582
#>  [697] 2001.591 2001.601 2001.610 2001.620 2001.629 2001.639 2001.648 2001.657
#>  [705] 2001.667 2001.676 2001.686 2001.695 2001.705 2001.714 2001.724 2001.733
#>  [713] 2001.743 2001.752 2001.762 2001.771 2001.780 2001.790 2001.799 2001.809
#>  [721] 2001.818 2001.828 2001.837 2001.847 2001.856 2001.866 2001.875 2001.884
#>  [729] 2001.894 2001.903 2001.913 2001.922 2001.932 2001.941 2001.951 2001.960
#>  [737] 2001.970 2001.979 2001.988 2001.998 2002.007 2002.017 2002.026 2002.036
#>  [745] 2002.045 2002.055 2002.064 2002.074 2002.083 2002.093 2002.102 2002.111
#>  [753] 2002.121 2002.130 2002.139 2002.149 2002.158 2002.168 2002.177 2002.186
#>  [761] 2002.196 2002.205 2002.215 2002.224 2002.233 2002.243 2002.252 2002.261
#>  [769] 2002.271 2002.280 2002.290 2002.299 2002.308 2002.318 2002.327 2002.336
#>  [777] 2002.346 2002.355 2002.365 2002.374 2002.383 2002.393 2002.402 2002.411
#>  [785] 2002.421 2002.430 2002.440 2002.449 2002.458 2002.468 2002.477 2002.485
#>  [793] 2002.494 2002.503 2002.512 2002.521 2002.530 2002.539 2002.549 2002.558
#>  [801] 2002.568 2002.577 2002.587 2002.596 2002.605 2002.615 2002.624 2002.634
#>  [809] 2002.643 2002.653 2002.662 2002.672 2002.681 2002.691 2002.700 2002.710
#>  [817] 2002.719 2002.729 2002.738 2002.748 2002.757 2002.767 2002.776 2002.786
#>  [825] 2002.795 2002.805 2002.814 2002.824 2002.833 2002.843 2002.852 2002.862
#>  [833] 2002.871 2002.880 2002.890 2002.899 2002.909 2002.918 2002.928 2002.937
#>  [841] 2002.947 2002.956 2002.965 2002.975 2002.984 2002.993 2003.003 2003.012
#>  [849] 2003.021 2003.031 2003.040 2003.049 2003.059 2003.068 2003.077 2003.086
#>  [857] 2003.096 2003.105 2003.114 2003.124 2003.133 2003.142 2003.152 2003.161
#>  [865] 2003.167 2003.173 2003.182 2003.191 2003.201 2003.210 2003.219 2003.229
#>  [873] 2003.238 2003.247 2003.257 2003.266 2003.275 2003.285 2003.294 2003.303
#>  [881] 2003.313 2003.322 2003.331 2003.341 2003.350 2003.359 2003.369 2003.378
#>  [889] 2003.388 2003.397 2003.406 2003.416 2003.425 2003.434 2003.444 2003.453
#>  [897] 2003.462 2003.472 2003.481 2003.490 2003.500 2003.509 2003.519 2003.528
#>  [905] 2003.538 2003.547 2003.557 2003.566 2003.575 2003.585 2003.594 2003.604
#>  [913] 2003.613 2003.623 2003.632 2003.642 2003.651 2003.661 2003.670 2003.679
#>  [921] 2003.689 2003.698 2003.708 2003.717 2003.726 2003.735 2003.745 2003.754
#>  [929] 2003.763 2003.772 2003.782 2003.791 2003.800 2003.809 2003.818 2003.828
#>  [937] 2003.837 2003.846 2003.855 2003.865 2003.874 2003.883 2003.892 2003.901
#>  [945] 2003.911 2003.920 2003.929 2003.938 2003.948 2003.957 2003.966 2003.976
#>  [953] 2003.985 2003.994 2004.003 2004.013 2004.022 2004.032 2004.041 2004.050
#>  [961] 2004.060 2004.069 2004.078 2004.088 2004.097 2004.106 2004.116 2004.125
#>  [969] 2004.135 2004.144 2004.154 2004.163 2004.173 2004.182 2004.192 2004.201
#>  [977] 2004.210 2004.220 2004.229 2004.239 2004.248 2004.256 2004.265 2004.274
#>  [985] 2004.283 2004.292 2004.301 2004.309 2004.317 2004.326 2004.334 2004.343
#>  [993] 2004.351 2004.360 2004.368 2004.377 2004.385 2004.393 2004.402 2004.410
#> [1001] 2004.419 2004.428 2004.437 2004.446 2004.455 2004.464 2004.473 2004.482
#> [1009] 2004.491 2004.500 2004.509 2004.519 2004.528 2004.538 2004.547 2004.556
#> [1017] 2004.566 2004.575 2004.585 2004.594 2004.603 2004.613 2004.622 2004.632
#> [1025] 2004.641 2004.650 2004.660 2004.669 2004.679 2004.688 2004.698 2004.707
#> [1033] 2004.716 2004.726 2004.735 2004.745 2004.754 2004.763 2004.773 2004.782
#> [1041] 2004.792 2004.801 2004.810 2004.820 2004.829 2004.839 2004.848 2004.857
#> [1049] 2004.867 2004.876 2004.884 2004.892 2004.900 2004.910 2004.919 2004.929
#> [1057] 2004.938 2004.947 2004.957 2004.966 2004.975 2004.985 2004.994 2005.004
#> [1065] 2005.013 2005.022 2005.032 2005.041 2005.051 2005.060 2005.067 2005.075
#> [1073] 2005.084 2005.092 2005.100 2005.108 2005.116 2005.125 2005.134 2005.143
#> [1081] 2005.152 2005.161 2005.170 2005.179 2005.188 2005.198 2005.207 2005.216
#> [1089] 2005.224 2005.232 2005.241 2005.249 2005.258 2005.266 2005.273 2005.280
#> [1097] 2005.287 2005.296 2005.306 2005.315 2005.324 2005.333 2005.343 2005.352
#> [1105] 2005.361 2005.370 2005.379 2005.387 2005.396 2005.405 2005.413 2005.422
#> [1113] 2005.431 2005.437 2005.443 2005.450 2005.457 2005.466 2005.475 2005.484
#> [1121] 2005.493 2005.500 2005.509 2005.519 2005.528 2005.538 2005.547 2005.556
#> [1129] 2005.566 2005.575 2005.584 2005.594 2005.603 2005.613 2005.622 2005.631
#> [1137] 2005.641 2005.650 2005.660 2005.669 2005.678 2005.688 2005.697 2005.707
#> [1145] 2005.716 2005.725 2005.735 2005.744 2005.753 2005.763 2005.772 2005.782
#> [1153] 2005.791 2005.800 2005.810 2005.819 2005.829 2005.838 2005.847 2005.856
#> [1161] 2005.866 2005.875 2005.884 2005.893 2005.902 2005.911 2005.920 2005.929
#> [1169] 2005.938 2005.947 2005.956 2005.965 2005.973 2005.983 2005.992 2006.002
#> [1177] 2006.011 2006.020 2006.030 2006.039 2006.048 2006.057 2006.066 2006.075
#> [1185] 2006.085 2006.094 2006.103 2006.112 2006.121 2006.130 2006.139 2006.149
#> [1193] 2006.158 2006.167 2006.176 2006.185 2006.194 2006.204 2006.213 2006.222
#> [1201] 2006.231 2006.241 2006.250 2006.259 2006.268 2006.278 2006.287 2006.296
#> [1209] 2006.305 2006.314 2006.324 2006.333 2006.342 2006.351 2006.360 2006.370
#> [1217] 2006.379 2006.388 2006.397 2006.406 2006.416 2006.425 2006.434 2006.443
#> [1225] 2006.451 2006.460 2006.469 2006.477 2006.486 2006.495 2006.500 2006.508
#> [1233] 2006.517 2006.525 2006.533 2006.542 2006.550 2006.559 2006.569 2006.578
#> [1241] 2006.587 2006.596 2006.606 2006.615 2006.624 2006.634 2006.643 2006.652
#> [1249] 2006.661 2006.671 2006.680 2006.689 2006.699 2006.708 2006.717 2006.727
#> [1257] 2006.736 2006.745 2006.754 2006.764 2006.773 2006.782 2006.792 2006.801
#> [1265] 2006.810 2006.820 2006.829 2006.838 2006.848 2006.857 2006.867 2006.876
#> [1273] 2006.886 2006.895 2006.904 2006.914 2006.923 2006.933 2006.942 2006.952
#> [1281] 2006.961 2006.971 2006.980 2006.989 2006.999 2007.008 2007.018 2007.027
#> [1289] 2007.037 2007.046 2007.056 2007.065 2007.075 2007.084 2007.093 2007.103
#> [1297] 2007.112 2007.122 2007.131 2007.141 2007.150 2007.160 2007.169 2007.179
#> [1305] 2007.188 2007.197 2007.207 2007.216 2007.226 2007.235 2007.245 2007.254
#> [1313] 2007.264 2007.273 2007.283 2007.292 2007.301 2007.311 2007.320 2007.330
#> [1321] 2007.339 2007.349 2007.358 2007.368 2007.377 2007.387 2007.396 2007.405
#> [1329] 2007.415 2007.424 2007.434 2007.443 2007.453 2007.462 2007.472 2007.481
#> [1337] 2007.491 2007.500 2007.509 2007.518 2007.528 2007.537 2007.546 2007.555
#> [1345] 2007.565 2007.574 2007.583 2007.592 2007.601 2007.611 2007.620 2007.629
#> [1353] 2007.638 2007.648 2007.657 2007.666 2007.675 2007.684 2007.694 2007.703
#> [1361] 2007.712 2007.721 2007.722 2007.732 2007.741 2007.751 2007.760 2007.770
#> [1369] 2007.779 2007.789 2007.798 2007.808 2007.817 2007.827 2007.836 2007.846
#> [1377] 2007.855 2007.865 2007.874 2007.884 2007.893 2007.903 2007.912 2007.922
#> [1385] 2007.931 2007.940 2007.950 2007.959 2007.969 2007.978 2007.988 2007.997
#> [1393] 2008.007 2008.016 2008.026 2008.035 2008.045 2008.054 2008.064 2008.073
#> [1401] 2008.083 2008.092 2008.102 2008.111 2008.121 2008.130 2008.140 2008.149
#> [1409] 2008.159 2008.168 2008.178 2008.187 2008.197 2008.206 2008.215 2008.225
#> [1417] 2008.234 2008.244 2008.253 2008.263 2008.272 2008.282 2008.291 2008.300
#> [1425] 2008.310 2008.319 2008.328 2008.337 2008.346 2008.356 2008.365 2008.374
#> [1433] 2008.383 2008.392 2008.402 2008.411 2008.420 2008.429 2008.438 2008.448
#> [1441] 2008.457 2008.466 2008.475 2008.483 2008.492 2008.500 2008.509 2008.519
#> [1449] 2008.528 2008.538 2008.547 2008.557 2008.566 2008.576 2008.585 2008.595
#> [1457] 2008.604 2008.613 2008.623 2008.632 2008.642 2008.651 2008.661 2008.670
#> [1465] 2008.680 2008.689 2008.698 2008.708 2008.717 2008.727 2008.736 2008.744
#> [1473] 2008.753 2008.763 2008.772 2008.782 2008.791 2008.801 2008.810 2008.820
#> [1481] 2008.829 2008.839 2008.848 2008.857 2008.867 2008.876 2008.886 2008.895
#> [1489] 2008.905 2008.914 2008.924 2008.933 2008.943 2008.952 2008.961 2008.971
#> [1497] 2008.980 2008.990 2008.999 2009.009 2009.017 2009.026 2009.034 2009.043
#> [1505] 2009.052 2009.060 2009.069 2009.077 2009.087 2009.096 2009.106 2009.115
#> [1513] 2009.124 2009.134 2009.143 2009.153 2009.162 2009.171 2009.181 2009.190
#> [1521] 2009.200 2009.209 2009.218 2009.228 2009.237 2009.246 2009.256 2009.265
#> [1529] 2009.275 2009.284 2009.293 2009.303 2009.312 2009.322 2009.331 2009.340
#> [1537] 2009.350 2009.359 2009.369 2009.378 2009.387 2009.397 2009.406 2009.415
#> [1545] 2009.425 2009.434 2009.444 2009.453 2009.462 2009.472 2009.481 2009.491
#> [1553] 2009.500 2009.507 2009.515 2009.522 2009.532 2009.541 2009.550 2009.560
#> [1561] 2009.569 2009.578 2009.588 2009.597 2009.606 2009.616 2009.625 2009.634
#> [1569] 2009.644 2009.653 2009.662 2009.672 2009.681 2009.690 2009.700 2009.709
#> [1577] 2009.718 2009.728 2009.737 2009.746 2009.756 2009.765 2009.774 2009.784
#> [1585] 2009.793 2009.802 2009.812 2009.821 2009.830 2009.840 2009.849 2009.858
#> [1593] 2009.868 2009.877 2009.886 2009.896 2009.905 2009.914 2009.924 2009.933
#> [1601] 2009.942 2009.952 2009.961 2009.970 2009.980 2009.989 2009.998 2010.008
#> [1609] 2010.017 2010.026 2010.036 2010.045 2010.055 2010.064 2010.073 2010.083
#> [1617] 2010.092 2010.101 2010.111 2010.120 2010.130 2010.139 2010.149 2010.158
#> [1625] 2010.168 2010.177 2010.187 2010.196 2010.206 2010.215 2010.225 2010.234
#> [1633] 2010.244 2010.253 2010.263 2010.272 2010.282 2010.291 2010.301 2010.310
#> [1641] 2010.320 2010.329 2010.339 2010.348 2010.358 2010.367 2010.377 2010.386
#> [1649] 2010.396 2010.405 2010.415 2010.424 2010.434 2010.443 2010.453 2010.462
#> [1657] 2010.472 2010.481 2010.491 2010.500 2010.509 2010.519 2010.528 2010.537
#> [1665] 2010.547 2010.556 2010.565 2010.575 2010.584 2010.593 2010.603 2010.612
#> [1673] 2010.621 2010.631 2010.640 2010.649 2010.659 2010.668 2010.677 2010.687
#> [1681] 2010.696 2010.705 2010.715 2010.724 2010.733 2010.743 2010.752 2010.761
#> [1689] 2010.771 2010.780 2010.789 2010.799 2010.808 2010.817 2010.827 2010.836
#> [1697] 2010.845 2010.855 2010.864 2010.871 2010.879 2010.888 2010.898 2010.907
#> [1705] 2010.917 2010.926 2010.936 2010.945 2010.955 2010.964 2010.973 2010.983
#> [1713] 2010.992 2011.002 2011.011 2011.021 2011.030 2011.040 2011.049 2011.058
#> [1721] 2011.068 2011.077 2011.087 2011.096 2011.106 2011.115 2011.124 2011.134
#> [1729] 2011.143 2011.153 2011.162 2011.172 2011.181 2011.191 2011.200 2011.209
#> [1737] 2011.219 2011.228 2011.238 2011.247 2011.256 2011.266 2011.275 2011.284
#> [1745] 2011.294 2011.303 2011.313 2011.322 2011.331 2011.341 2011.350 2011.359
#> [1753] 2011.369 2011.378 2011.388 2011.397 2011.406 2011.416 2011.425 2011.434
#> [1761] 2011.444 2011.453 2011.463 2011.472 2011.481 2011.491 2011.500 2011.509
#> [1769] 2011.519 2011.528 2011.538 2011.547 2011.557 2011.566 2011.576 2011.585
#> [1777] 2011.595 2011.604 2011.614 2011.623 2011.633 2011.642 2011.652 2011.661
#> [1785] 2011.671 2011.680 2011.690 2011.699 2011.709 2011.718 2011.728 2011.737
#> [1793] 2011.747 2011.756 2011.766 2011.775 2011.785 2011.794 2011.804 2011.813
#> [1801] 2011.823 2011.832 2011.842 2011.851 2011.861 2011.870 2011.880 2011.889
#> [1809] 2011.899 2011.908 2011.918 2011.927 2011.937 2011.946 2011.956 2011.965
#> [1817] 2011.975 2011.984 2011.994 2012.003 2012.013 2012.022 2012.032 2012.041
#> [1825] 2012.051 2012.060 2012.070 2012.079 2012.089 2012.098 2012.108 2012.117
#> [1833] 2012.127 2012.136 2012.146 2012.155 2012.165 2012.174 2012.184 2012.193
#> [1841] 2012.203 2012.212 2012.222 2012.231 2012.241 2012.250 2012.260 2012.269
#> [1849] 2012.279 2012.288 2012.298 2012.307 2012.317 2012.326 2012.336 2012.345
#> [1857] 2012.355 2012.364 2012.374 2012.383 2012.393 2012.402 2012.412 2012.421
#> [1865] 2012.431 2012.440 2012.450 2012.459 2012.469 2012.478 2012.488 2012.497
#> [1873] 2012.507 2012.516 2012.526 2012.535 2012.545 2012.554 2012.564 2012.573
#> [1881] 2012.583 2012.592 2012.602 2012.611 2012.620 2012.630 2012.639 2012.648
#> [1889] 2012.658 2012.667 2012.677 2012.686 2012.695 2012.705 2012.714 2012.723
#> [1897] 2012.733 2012.742 2012.752 2012.761 2012.770 2012.780 2012.789 2012.798
#> [1905] 2012.808 2012.817 2012.826 2012.836 2012.845 2012.855 2012.864 2012.873
#> [1913] 2012.883 2012.892 2012.901 2012.911 2012.920 2012.930 2012.939 2012.948
#> [1921] 2012.958 2012.967 2012.976 2012.986 2012.995 2013.004 2013.014 2013.023
#> [1929] 2013.032 2013.042 2013.051 2013.061 2013.070 2013.079 2013.089 2013.098
#> [1937] 2013.107 2013.117 2013.126 2013.135 2013.145 2013.154 2013.163 2013.173
#> [1945] 2013.182 2013.191 2013.201 2013.210 2013.219 2013.229 2013.238 2013.248
#> [1953] 2013.257 2013.266 2013.276 2013.285 2013.294 2013.304 2013.313 2013.322
#> [1961] 2013.332 2013.341 2013.350 2013.360 2013.369 2013.378 2013.388 2013.397
#> [1969] 2013.406 2013.416 2013.425 2013.435 2013.444 2013.453 2013.463 2013.472
#> [1977] 2013.481 2013.491
pdf("Figures/s5_90.pdf",6,6)
plot(s5a)
plot(s5b)
plot(s5c)
dev.off()
#> png 
#>   2
s5a

s5b

s5c

Hyppeparameter pairs

s6 <- plot_hyperpar_pairs(out)
pdf("Figures/s6_90.pdf",6,6)
plot(s6)
#> Warning: The dot-dot notation (`..density..`) was deprecated in ggplot2 3.4.0.
#> ℹ Please use `after_stat(density)` instead.
#> ℹ The deprecated feature was likely used in the ResistPhy package.
#>   Please report the issue to the authors.
dev.off()
#> png 
#>   2
s6

s7a <- plot_epi_dynamic(out,1)
s7b <- plot_rt(out,1)
pdf("Figures/s7_90.pdf",6,6)
plot(s7a)
plot(s7b)
dev.off()
#> png 
#>   2
s7a

s7b